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Abstract 

Intense,  ultrashort  laser  pulses  propagating  in  the  atmosphere  have  been  observed 
to  emit  sub-THz  electromagnetic  pulses  (EMP).  The  purpose  of  this  paper  is  to  analyze 
EMP  generation  from  the  interaction  of  ultrashort  laser  pulses  with  air  and  with  dielectric 
surfaces  and  to  determine  the  efficiency  of  conversion  of  laser  energy  to  EMP  energy.  In 
our  self-consistent  model  the  laser  pulse  partially  ionizes  the  medium,  forms  a  plasma 
filament,  and  through  the  ponderomotive  forces  associated  with  the  laser  pulse,  drives 
plasma  currents  which  are  the  source  of  the  EMP.  The  propagating  laser  pulse  evolves 
under  the  influence  of  diffraction,  Kerr  focusing,  plasma  defocusing,  and  energy 
depletion  due  to  electron  collisions  and  ionization.  Collective  effects  and  recombination 
processes  are  also  included  in  the  model.  The  duration  of  the  EMP  in  air,  at  a  fixed  point, 
is  found  to  be  a  few  hundred  femtoseconds,  i.e.,  on  the  order  of  the  laser  pulse  duration 
plus  the  electron  collision  time.  For  steady  state  laser  pulse  propagation  the  flux  of  EMP 
energy  is  non-radiative  and  axially  directed.  Radiative  EMP  energy  is  present  only  for 
non-steady  state  or  transient  laser  pulse  propagation.  The  analysis  also  considers  the 
generation  of  EMP  on  the  surface  of  a  dielectric  on  which  an  ultrashort  laser  pulse  is 
incident.  For  typical  laser  parameters,  the  power  and  energy  conversion  efficiency  from 
laser  radiation  to  EMP  radiation  in  both  air  and  from  dielectric  surfaces  is  found  to  be 
extremely  small,  <  10‘8.  Results  of  full  scale,  self-consistent,  numerical  simulations  of 
atmospheric  and  dielectric  surface  EMP  generation  are  presented.  A  recent  experiment 
on  atmospheric  EMP  generation  is  also  simulated. 
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I.  Introduction 

Propagation  of  intense,  ultrashort  laser  pulses  in  the  atmosphere  has  demonstrated  a  wide 
range  of  phenomena  that  may  have  a  variety  of  applications.  The  observations  include 
formation  and  stable,  long-distance  propagation  of  optical  and  plasma  filaments,  ‘white’ 
light  generation  and  emission  of  sub-THz  electromagnetic  pulses  (EMP)  [1-18]. 

Terahertz  radiation  has  a  wide  range  of  applications  in  areas  such  as  chemical 
identification,  biomedical  diagnostics  and  dielectric  measurements.  In  addition, 
tomographic  imaging  using  THz  radiation  has  potential  applications  in  biomedical 
research,  security  screening  and  remote  sensing.  Terahertz  radiation  is  also  generated  by 
backward  wave  oscillators  [19],  gyrotrons  [20],  ffee-electron  lasers  [21],  difference- 
frequency  mixing  in  nonlinear  media  [22],  photoconductive  switches  [23],  tunable  optical 
parametric  oscillators[24],  semiconductor  lasers  [25],  optically-pumped  THz  lasers  [26], 
DC  to  AC  converters  [27]  and  Cherenkov  wakes  in  magnetized  plasmas  [28]. 

Recently,  it  has  been  observed  that  a  broadband  electromagnetic  pulse  (EMP)  of 
terahertz  radiation  can  be  generated  within  the  plasma  filaments  formed  by  an  ultra-short 
(~  100  fsec),  high-intensity  (~10I3-1014  W/cm2)  laser  pulse  propagating  in  air  [14,  17],  In 
these  experiments  the  laser  pulse  was  sufficiently  intense  to  create  a  plasma  filament 
through  multiphoton  ionization  of  air.  The  plasma  density  within  the  filament  was 
observed  to  be  ~  1016  cm  3 .  This  density  corresponds  to  a  plasma  wavelength  of  ~ 

330 /urn  and  a  plasma  frequency  of  ~  0.9  THz. 

In  this  paper  we  analyze  and  numerically  simulate  the  EMP  produced  by  the 
interaction  of  ultrashort  laser  pulses  with  air  and  dielectric  surfaces,  and  determine  the 
conversion  efficiency  from  the  laser  pulse  energy  to  EMP.  Atmospheric  EMP  is 
generated  when  a  short,  intense  laser  pulse  ionizes  the  air  and  forms  a  plasma  filament. 
The  plasma  filament  extends  behind  and  its  front  travels  with  the  laser  pulse  as  indicated 
in  Fig.  1 .  The  duration  of  the  plasma  current  which  drives  the  EMP  is  limited  by  the 
electron  collision  frequency.  For  a  plasma  density  of  ne  ~  1016  cm-3 ,  the  electron 

collision  frequency  is  ve  ~  5xl012  sec-1  (  for  Te  =  leV)  and  the  duration  of  the  plasma 
current  is  ~  1  /  ve  ~  200  fsec .  The  duration  of  the  EMP  in  air,  at  a  fixed  point,  is  a  few 
hundred  femtoseconds,  i.e.,  the  laser  pulse  duration  plus  the  electron  collision  time.  The 
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spatial  gradients  associated  with  the  laser  pulse  envelope  as  well  as  electron  collisions 
lead  to  ponderomotive  forces  on  the  plasma  electrons.  These  forces  cause  the  electrons 
to  oscillate,  setting  up  plasma  currents.  The  axial  and  transverse  components  of  plasma 
current  pulse  have  velocities  close  to  the  group  velocity  of  the  laser  pulse.  This  is  similar 
to  the  laser  driven  excitation  of  plasma  waves  in  the  laser  wakefield  accelerator  (LWFA) 
[29].  In  the  LWFA  plasma  oscillations  occur  at  the  local  plasma  frequency.  In  general 
the  various  field  components  of  the  plasma  oscillations  are  coupled,  as  are  the  axial  and 
transverse  plasma  currents.  These  induced  plasma  currents,  which  travel  with  the  laser 
pulse,  are  the  source  of  the  atmospheric  EMP.  The  axial  extent  of  the  EMP  is  limited, 
and  consequently,  the  spectrum  of  the  EMP  pulse  can  be  broad.  However,  the 
characteristic  frequency  of  the  EMP  measured  outside  of  the  plasma  filament  is 
determined  by  the  laser  pulse  duration  and  not  the  local  plasma  frequency. 

The  results  of  our  analysis  for  atmospheric  EMP  generation  can  be  summarized  as 
follows.  For  steady  state  propagation,  in  which  the  laser  pulse  and  plasma  filament  are 
static  in  the  group  velocity  frame  of  the  laser  pulse,  the  plasma  current  source  travels  with 
a  velocity  equal  to  the  laser  pulse  group  velocity.  Since  the  group  velocity  is  necessarily 
less  than  the  speed  of  light,  there  is  no  EMP  radiation  generated.  The  absence  of  radial 
EMP  emission  in  the  steady  state  is  one  of  several  fundamental  differences  between  our 
results  and  those  of  Refs.  [30]  and  [31],  see  also  [32].  The  absence  of  EMP  in  the  steady 
state  is  analogous  to  the  lack  of  Cherenkov  radiation  when  the  velocity  of  a  charged 
particle  is  less  than  the  velocity  of  light  in  the  medium.  The  radial  component  of  the 
plasma  current  in  the  filament  can  lead  to  surface  charges  on  the  plasma-air  boundary. 

The  surface  charge  density  is  strongly  dependent  on  the  collisionality  of  the  plasma.  In 
the  weakly  collisional  regime  the  surface  charge  leads  to  greatly  reduced  EMP  fields 
beyond  the  filament. 

Under  conditions  of  non-steady  state  (transient)  propagation,  the  plasma  currents 
can  generate  EMP  having  a  radial  energy  flux.  For  non-steady  state  pulse  propagation, 
the  plasma  currents  can  have  Fourier  components  with  phase  velocities  greater  than  the 
laser  pulse  group  velocity.  These  Fourier  components  can  generate  radiative  EMP, 
analogous  to  Cherenkov  radiation.  A  condition  for  generating  radiative  EMP  in  the 
transient  regime  is  derived  and  the  differing  characteristics  of  EMP  in  the  steady  state  and 
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transient  propagation  regimes  are  discussed  in  detail  and  illustrated  with  simulations.  For 
typical  atmospheric  conditions  and  laser  parameters  it  is  found  that  the  efficiency  of  the 
EMP  generation  from  laser-atmospheric  interaction,  i.e.,  the  ratio  of  the  EMP  energy  to 
the  laser  pulse  energy,  is  extremely  small. 

In  addition  to  atmospheric  EMP  generation  we  consider  EMP  from  the  interaction 
of  short  laser  pulses  with  dielectric  surfaces.  The  source  of  the  EMP  on  the  surfaces  of 
the  dielectric  is  the  oscillating  plasma  sheet  formed  by  the  ionization  process  as  well  as 
an  optical  rectification  process.  Our  analysis  of  EMP  generation  on  the  surfaces  of  a 
dielectric  indicates  that  the  conversion  efficiency  is  extremely  small. 

This  paper  is  organized  as  follows.  The  general  equations  describing  the 
evolution  of  the  driving  laser  pulse  and  EMP  in  air  are  derived  in  Sec.  II.  Also  in  Sec.  II, 
an  estimate  of  the  EMP  power  conversion  efficiency  is  derived  and  an  analysis  of  the 
steady  state  and  transient  propagation  regimes  presented.  In  Sec.  Ill  EMP  generation 
from  a  dielectric  surface  on  which  a  laser  pulse  is  incident  is  analyzed.  Section  IV 
presents  self-consistent  numerical  simulations  showing  laser  pulse  propagation,  plasma 
filament  formation,  and  EMP  generation  in  air  and  on  dielectric  surfaces.  A  summary  of 
the  results  and  concluding  remarks  are  found  in  Sec.  V. 

II.  Electromagnetic  Pulse  Generation  in  Air 

In  this  section  we  derive  the  general  equations  for  the  generation  of  EMP  radiation 
from  laser-induced  plasma  currents  in  air.  For  the  parameters  of  interest,  the  driving  laser 
pulse  intensity  is  sufficiently  high  that  it  partially  ionizes  the  air  and  forms  a  plasma 
filament  behind  the  laser  pulse.  An  equation  for  the  amplitude  of  the  driving  laser  pulse 
is  derived  that  includes  the  effects  of  diffraction,  Kerr  focusing,  plasma  defocusing, 
electron  collisions  and  laser  pulse  energy  depletion  caused  by  ionization.  The  plasma 
currents  within  the  filament  are  driven  by  the  ponderomotive  forces  associated  with  the 
laser  pulse  and  are  described  by  a  set  of  cold  fluid  equations  which  include  the  effects  of 
multi-photon  ionization,  electron-neutral  elastic  collisions  and  electron 
recombination/attachment.  In  our  model  an  axially  symmetric  laser-plasma  configuration 
is  assumed  and  stimulated  Raman  scattering  of  the  laser  pulse  and  other  instabilities  are 
not  considered  [33,  34]. 
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a)  Driving  Laser  Pulse 

For  a  self-consistent  description  of  EMP  generation  the  spatial  and  temporal 
evolution  of  the  driving  laser  pulse  is  needed.  The  ponderomotive  forces,  which  are 
proportional  to  the  energy  density  of  the  driving  laser  pulse,  induce  low  frequency  plasma 
currents  which  in  turn  generate  EMP.  The  driving  laser  pulse  is  assumed  to  be  many 
optical  wavelengths  long  and  to  have  a  spot  size  that  is  large  compared  to  the  wavelength. 
The  transverse  laser  fields  are  given  by 


Ei(r,/)  = 

=  EL (r, t) exp(  i(k0  z  -  co0  t))ex  / 2+  c.c.. 

(la) 

Bi(r  ,0  = 

-  BL(r,t)exp(i(k0  z  -  0)ot))ey/2  +  c.c., 

(lb) 

where  a>0  is  the  carrier  frequency,  ko  =  (0o!c  is  the  wavenumber,  ex  (e^jis  a  unit 

vector  in  the  x  (y)  direction  and  c.c.  denotes  the  complex  conjugate.  In  Eq.(la)  EL  is  the 

slowly  varying  complex  electric  field  amplitude  satisfying  jiff1  d  EL  /  dt\  «  co0  and 

\e~lx  dELl dzj  «  k0 .  The  transverse  magnetic  field  amplitude  is 

Bl  =  [1  -  i(c/co0)(d/dz  +  c~'d/dt)]EL .  The  axial  electric  and  magnetic  fields  of 

magnitude  ~(c/ co0)\d  EL  /  dx\  and  ~  {c  I  coa)\d  EL/  d  y\,  respectively,  that  are  associated 

with  the  transverse  laser  fields  in  Eqs.(l),  are  small  and  are  neglected. 

Making  use  of  the  electron  cold  fluid  equations,  the  transverse  current  density 
J L  -  qne  v L  induced  by  the  laser  pulse  is  given  by 


dJL 

dt 


+  VeJL 


(2) 


where  ve  =  ven  +  vei  is  the  electron  collision  frequency,  y  is  the  electron-neural 
collision  frequency,  vei  is  the  electron-ion  collision  frequency, 

a>p(r,t)  =  (Aitq2  ne(r,t)/ m)l/2  is  the  local  plasma  frequency,  q  is  the  electronic  charge, 


\L  is  the  quiver  velocity  in  the  laser  field  and  ne  is  the  electron  density  in  the  plasma 
filament  which  is  formed  by  photoionization.  It  is  assumed  that  the  ions  are  immobile  on 
the  time  scale  of  interest.  In  general,  the  electron  density  ne  includes  contributions  from 
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plasma  oscillations  as  well  as  the  ionization  process.  The  solution  of  Eq.(2)  for  the  first 
order  transverse  current  density  induced  by  the  laser  field  is 

Jr(r»0  =  JL(r,t)exp(i(k0z  -a0t))ex/2  +  c.c.,  (3) 

where  J  L  =  (it  Anco0 )  [  1  -  (i/  a>0)(d/dt  +  ve)]  ( co2pEL  ) ,  where  we  have  assumed  that 

The  electron-neutral  collision  frequency  is  given  by  ven  =  nn  er0  (kB  Te/m)U2 
where  nn  -  nno  -  ne  is  the  neutral  density  (singly  ionized  plasma),  nn0  is  the  ambient 
neutral  density,  kB  is  the  Boltzmann  constant  and  a0  «  5  x  10“15  cm2  is  a  typical  value 
for  the  scattering  cross  section  [35].  In  practical  units  the  expression  can  be  written  as 
ve«[sec  ']  ~  2xl0'7  n(I[cm"3](re[eV])1/2 .  For  a  neutral  air  density  of 

n„  =  n„0  =  2.7xl019cm  3  and  temperature  Te  =  1  eV ,  the  electron  collision  frequency 
is  vm  «  5.7xl012  sec-1 .  Note  that  for  a  plasma  density  of  ne  ~  1016  cm-3 ,  the  electron- 
neutral  collision  frequency  and  plasma  frequency  are  comparable,  ven  ~  a>p .  Taking  the 
Coulomb  logarithm  to  be  equal  to  10,  the  electron-ion  collision  frequency  [35]  is 
ve(.[sec  ‘]«3xl0  5  n-t  [cm-3] (^[eV])-372 ,  for  singly  ionized  plasma,  i.e.,  ne  =  ni .  The 
electron-neural  and  electron-ion  collision  frequency  become  equal  when  the  electron 
density  reaches  the  value  ne  *  1.8xl017  T2[eV]  for  singly  ionized  plasma. 
i)  Plasma  Current 

The  propagation  of  the  driving  laser  pulse  is  strongly  influence  by  the  plasma 
filament.  The  plasma  dynamics  are  described  by  the  non-relativistic,  cold  fluid  equations 

dn 

~zf  +  V  •  (ne  v)  =  vion  n„  -  precom  n]  -flamch  ne  n] ,  (4a) 


— +  (v-V)v  =  —  E  + 
dt  m  . 


vxB 


-  V.  V  - 


+  V-(«ev)  \/ne. 


where,  v  is  the  electron  fluid  velocity,  vitm  is  the  ionization  rate,  firecom  is  the 
recombination  coefficient,  (5attach  is  the  attachment  coefficient,  nu  =  nno  -  ne , 
ve  =  Ve„  +  v ei  an^  the  ions  are  assumed  to  be  stationary.  In  obtaining  Eq.(4)  we  have 
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assumed  that  the  electrons  immediately  after  ionization  have  zero  velocity,  i.e.,  the 
electrons  are  bom  at  rest  and  above  threshold  ionization  effects  are  ignored. 

For  the  parameters  of  interest  here  the  ionization  process  is  principally  due  to 
multi-photon  ionization  and  the  corresponding  rate  can  be  fitted  analytically  by  [36] 


where  the  integer  l  <  Int(Ui0H  /hco0  +  1)  refers  to  the  effective  number  of  photons 


absorbed  by  the  molecule  in  the  ionization  process,  U ion  is  the  ionization  energy,  IL  is 


the  laser  intensity  and  Iup  -tio)2o  I  oMP ,  with  aMP  an  empirically  determined  ionization 


cross  section.  At  laser  intensities  of  ~1014  W/cm2or  above  tunneling  ionization  becomes 
important.  For  a  laser  wavelength  of  XQ  =  0.8  pm  (ho)0  =1.54eV)  in  air,  the  multi¬ 
photon  ionization  order  of  02  is  l  =  8  ,  and  the  characteristic  intensity  is  determined 
empirically  to  be  lMP  ~  1014  W/cm2  for  pulses  with  durations  ~  100  fsec  [16].  Since 
the  intensity  ratio  in  Eq.(5)  is  raised  to  the  eighth  power,  the  ionization  rate  is  an 
exceedingly  sensitive  function  of  the  laser  intensity  and  the  value  of  the  parameter  IMP . 
The  scaling  of  the  ionization  rate  in  Eq.(5)  with  laser  intensity  is  the  same  as  that  given  in 
Ref.  [37].  Since  the  numerical  value  of  the  Keldysh  formula  in  [37]  differs  from  the 
experimental  value,  we  use  the  empirical  cross-section  quoted  above.  A  similar  approach 
was  adopted  in  Ref.  [38]  to  obtain  an  empirical  formula  for  the  multi -photon  ionization 
rate.  For  a  laser  intensity  of  IL  =  7  x  1013  W/cm2 ,  laser  wavelength  of 
XQ  =0.8  fim  ( co0  =  2.4xl015  sec-1)  and  ionization  order  of  t  =  8,  the  ionization  rate  is 
vion  ~  2x1 0U  sec-1 .  The  ionization  rate,  however,  is  an  extremely  sensitive  function  of 


the  laser  intensity.  The  recombination  time  is  given  by  trecom  =  1  /  Precom  ne .  Taking 
Aw™  [cm3 /sec]  *  2xl0-8  (TJeV])-172  [14],  ne  =  1016  cm-3 ,  and  Te  =leV,the 


recombination  time  Trecom  ~  5  nsec  The  attachment  coefficient  Patlach  of  electrons  on 
oxygen  molecules  is  discussed  in  Ref  [39].  For  an  electron  temperature  of  ~  1  eV  the 
attachment  coefficient  is  paltach  ~  3.4x1 0-31  cm6/sec  and  the  attachment  time  is 
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T attach  =  1 !  P attach  nl  ~  4  nsec.  The  duration  of  the  plasma  filament  is  on  the  order  of  the 
shorter  of  the  two  times  scales,  zrecom  and  ratlach,  i.e,  rmn[r recom,  r  attach] . 

Combining  Eqs.(4a)  and  (4b)  results  in  an  equation  for  the  electron  plasma 
current  density,  J  =  qne\  ,  given  by  [40] 


SJ 

dt 


00  n 

—  E  +  S 
An 


(6) 


where  S  (r,z,t)  is  the  nonlinear  driving  term.  The  first  term  on  the  right-hand-side  of 
Eq.(6)  is  associated  with  collective  effects  while  the  second  term  represents  the  driving 
ponderomotive  force  due  to  the  laser  pulse.  The  duration  of  the  plasma  current  density  is 
limited  to  ~l/vg  ~  200fsecfor  ne  =10I6cm-3.  To  second  order  in  the  laser  field 
amplitude, 


S  =  -  <7~'(Ji-V)(Ji/«e).  (7) 

When  the  laser  spot  size  is  large  compared  to  the  plasma  filament  spot  size,  i.e., 
R0  »  rp ,  the  ponderomotive  source  term  in  Eq.(7)  can  be  simplified.  This  is  a  valid 

approximation  in  the  case  of  higher  order  multi -photon  ionization  where  the  ionization 
rate  is  a  strong  function  of  laser  intensity.  In  this  limit  the  axial  component  of  the  source 
term  S2  is  given  by 


Sz(r,z,t)  w  - 


®  p  (c  ^5  0 


EL(r,z,t)  |2.  (8) 


Note  that  the  axial  ponderomotive  force  consists  of  a  contribution  from  the  spatial 
gradient  of  the  laser  pulse  envelope  as  well  as  contributions  from  electron  collisions  and 
changes  in  the  electron  density  (ionization).  The  transverse  component  of  the  source 
term,  for  the  laser  pulse  linearly  polarized  in  the  x  direction,  is 


Si 


q  co2(r,z,t)(  Q  „  2  dco2 

- - -  0  -  g _ lL 

16 nm  co]  '  dr  *  a2p  dx 


(9) 


where  e(.  is  a  unit  vector  in  the  radial  direction.  Note  that  the  transverse  part  of  the 
source  term  has  a  component  in  the  x  direction.  For  a  circularly  polarized  laser  pulse  the 


9 


axial  component  of  the  source  term  is  given  by  Eq.(8)  with  EL  replaced  with  4lEL  and 
the  transverse  component  is  purely  radial  and  given  by 


q  a](r,z,t)(  d  1  da] 


EL(r,z,t)\  e,. 


The  axially  symmetric  circularly  polarized  laser  field  has  been  assumed  to  be  of  the  form 
E£  -  El  exp  (i(k0z  -  a0t))(ex  -  iey)/2  +  c.c. 


ii)  Laser  Pulse  Amplitude 

To  obtain  the  ponderomotive  source  term  self  consistently,  an  equation  describing 
the  evolution  of  the  driving  laser  pulse  is  needed.  For  this  purpose  it  is  useful  to 
transform  to  the  group  velocity  frame  in  which  the  independent  variables  are  z  and 
t  =  t  -  z/vg  where  vg  denotes  the  linear  group  velocity  of  the  laser  pulse,  and  is 
defined  later.  In  terms  of  the  new  independent  variables  ( z ,  z )  the  derivatives  are  given 
by  d/dt  — »  d/dr  and  d/dz  — »  d/dz  -v~ld/dr.  In  the  variables  (z,r,r)  the  laser 

pulse  evolves  slowly,  i.e.,  is  nearly  stationary.  The  wave  equation  for  the  complex  laser 
amplitude  EL(r,z,r)  in  the  new  independent  variables  takes  the  form  [36] 


d££p-,z,r) 

dz 


__,2  col  nK 

V2  +  «  K 

-L  A 


4  7tc  n 


±-bl' 


&L(r,z,T) 


The  first  term  on  the  right  hand  side  of  Eq.(l  1)  is  the  Laplacian  operator  which  is 
responsible  for  diffractive  spreading  of  the  laser  beam.  The  second  term  is  proportional 
to  the  nonlinear  Kerr  refractive  index  nK  and  can  lead  to  nonlinear  self- focusing  [40,  41] 

The  terms  proportional  to  co2po  represent  the  effect  of  plasma  defocusing  and  collisions 

[36,  40].  The  last  term  describes  laser  pulse  energy  depletion  due  to  ionization  of  the  air 
molecules  with  effective  ionization  energy  Uion .  In  general  the  refractive  index  for 

electromagnetic  waves  in  a  medium  consists  of  contributions  from  the  linear  and 
nonlinear  polarization  field.  In  the  present  analysis  it  is  assumed  that  the  contribution  to 
the  nonlinear  refractive  index  is  proportional  to  the  laser  intensity,  i.e.,  Kerr  effect.  This 
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is  a  valid  approximation  for  in  a  gas  as  well  as  for  an  amorphous  dielectric.  For  the  laser 
pulse  this  nonlinear  contribution  to  the  index  leads  to  the  phenomena  of  nonlinear  Kerr 
self  focusing. 

Hi)  Source  Dependent  Expansion  Analysis  of  Laser  Pulse 

Equation  (11)  for  the  laser  pulse  amplitude  can  be  approximately  solved  by 
employing  the  source  dependent  expansion  (SDE)  method  [42].  In  this  method  the  laser 
field  is  assumed  to  have  a  Gaussian  transverse  profile 

\EL{r,ztr)  |  =  |£i(0,z,r)|exp(-r2/i?2(z,r)),  (12) 

where  R  is  the  laser  spot  size.  For  the  Gaussian  transverse  profile  the  laser  power  is 
Pl(z,t)  =  c\El{Q,z,t)\2  R\z,z)l\6.  (13) 


Employing  the  SDE  method  results  in  coupled  equations  for  the  spot  size  and  power, 
d2R(z,r)  4c2  1 


dz2 


dPL(z,T) 

dz 


R  (ZA) 


l_^(£lO  +  2££r.B.(0,,,r)  ,  ' 

PK  (t  +  1)! 


1  Ve  D  <  \  TT  2 

i  + 1  T z’ _  27  u‘°" Vion  n" R  ^ 5 


W{z,t),  (14a) 

(14b) 


where  PK  =  A]  !{2nnK  nn  / nno)  is  the  nonlinear  Kerr  focusing  power,  r  =  q2  / me 2  is 
the  classical  electron  radius  and 


W(z,t) 


R 


dz 


(r2u)  +  r2u2  ), 


U(z,t) 


(*z _  1  <  v. 

2 1  PL  dz  2^(^+l)2  co20  c  ’ 


dnen 

eo 

dz 


=  v.  n 

ton  n 


(15a) 

(15b) 

(15c) 


The  decrease  in  the  total  laser  pulse  power  as  a  function  of  propagation  distance  is 
due  to  both  collisional  absorption  and  ionization.  The  evolution  of  the  total  laser  pulse 
power  given  by  Eq.(14b)  can  be  rewritten  as 


dPiM 

dz 


V  R  coll 


+ 


PLM, 


(16) 


-'ion  ) 


where 
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(17a) 


-'coll 


(t+ 1) 


C0„ 


®po  Ve 


is  the  absorption  length  due  to  collisional  effects  and 

21  PL 


L  = 


”  U ion  Vion  "nR 


(17b) 


is  the  absorption  length  due  to  ionization  effects.  It  should  be  noted  that  due  to  the 
nonlinear  character  of  Eq.(16),  the  expressions  in  Eqs.(17)  are  at  best  local  estimates  for 
the  scale  lengths.  As  an  example,  consider  the  case  of  a  laser  pulse  having  a  power 
PL  =  0.1  TW ,  wavelength  Xo  =  0.8 pm,  pulse  duration  rL  =  lOOfsec,  intensity 


I L  =  5xl013  W/cm2  and  spot  size  R  ~  360 jum ,  propagating  in  air  with  a  neutral  density 
nno  =  2.7  xlO19  cm-3 ,  effective  ionization  energy  of  Uion  =  8eV  (i  =  8)  and  electron 
collision  frequency  of  ve  =  ven  +  vei  «  6xl012  sec'1 .  In  this  case  the  ionization  rate  is 
vion  ~  1010  sec"'  (IMP  =  1014W/cm2)  and  the  electron  density  is  ne  ~3xl016cm~3, 
{co0  / a>p0  ~  230) .  For  these  parameters  the  collisional  absorption  length  is  Lcoll  ~  25m 
while  the  ionization  absorption  length  is  Lion  ~  10  m  .  In  this  example,  the  laser  pulse 

propagation  distance  is  limited  by  ionization  absorption.  It  should  be  stressed  that  these 
characteristic  scale  lengths  are  very  sensitive  functions  of  the  laser  and  plasma 
parameters. 

iv)  Pulse  Group  Velocity 

The  group  velocity  is  given  by  \ g  =  c{n  +  co0  dn / da>0)~ 1 ,  where  n  is  the  total 
refractive  index  and  expressible  as 


no  ~ 


fa 


2 

po 


2c2 


2  col 


col  Rl 


+  nKlL- 


(18) 


i  1 2 

Here  R0  is  the  initial  spot  size  of  the  laser  pulse,  and  I L  =  cn0  \EL\  /8zr  is  the  time- 
averaged  laser  intensity.  The  filling  factor  /  is  the  ratio  of  the  plasma  filament  cross 
sectional  area  to  the  laser  pulse  cross  sectional  area.  Using  the  expression  for  the  plasma 
density  due  to  the  multi -photon  ionization  process  in  Eq.(5)  we  find  that  /  ~  (f +1)~‘ . 
The  filling  factor  is  necessary  for  calculating  the  correct  group  velocity  of  the  laser  pulse 
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in  the  presence  of  the  plasma  filament.  The  third  term  on  the  right  hand  side  of  Eq.(18) 
denotes  the  correction  to  the  refractive  index  due  to  the  finite  spot  size  of  the  laser  pulse 
and  the  fourth  term  represents  the  nonlinear  index  due  to  the  Kerr  effect.  Making  use  of 
the  expression  for  the  refractive  index,  the  group  velocity 
vg  =c/3g  =  \l(n+d(a>0  n)l d co0) is  expressible  as 


P* 


1  fa> 


2 

po 


2  c2 


no 


colR] 


nK  I L  • 


(19) 


In  deriving  Eq.(19),  it  is  assumed  that  the  variation  of  the  linear  index  and  Kerr  index 
with  frequency  is  negligible.  For  I L  =  1014  W/cm2,  R0  =  260 pm,  T0  =  0.8pm, 

ne  =  1016  cm  3,  nK  =  6xl0“19  cm2/W ,  and  l  =  8,  the  magnitudes  of  the  plasma,  finite 
spot  size,  and  optical  Kerr  contributions  to  J3g  are,  ~  3  x  10~7 ,  3  x  10"6 ,  and  ~  6  x  10'5 , 
respectively. 

b)  EMP  Equations  in  Air 

The  wave  equation  for  the  EMP  electric  field  E  in  air  which  is  driven  by  plasma 
currents  is  given  by 

n2  d2 E  4;r  <3  J 


VxVxE  + 


c2  dt1 


c2  dt  ’ 


(20) 


where  n0  is  the  linear  refractive  index  and  the  plasma  current  is  given  by  Eq.(6).  We  use 
a  cylindrical  coordinate  system  where  r,  0,  and  z  denote  the  radial  coordinate, 
azimuthal  angle,  and  axial  coordinate,  respectively,  and  assume  axial  symmetry,  i.e., 
d  /  dO  =  0  .  The  generation  of  EMP  from  the  surface  of  dielectrics  is  analyzed  in  Sec.  III. 

In  the  laser  pulse  group  velocity  frame,  the  axial  and  radial  components  of  the  EMP 
electric  field  are  given  by 


ld_ 
r  dr 

1 


f 


\ 


2  h 


V  drJ 
d2 


Krldr 2 


V  8  '  8 


1  d 
c 2  dr2 


2  a2 

vg  dzdr 


E  + 


+ 


1  d 


dz2 


KV8  dT 


Er  + 


_d_ 

dz 


1  9  /  >  4  n  dJ 

yEr  )  = 


r  dr 


f 


1  5 


vvg  dT 


d_ 

dz 


dE, 


c  dr 
4  k  dJ„ 


dr 


dr 


(21a) 


(21b) 


where  yg  =  (1  -  «2v2  / c2)'l/2 .  Using  Eq.(19)  y  is  found  to  be  given  by 
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rg  =  <t>o/(fa>2po  +  4 c2  !R2o  +  2a>lnKILy\  (22) 

For  X0  =  0.8 jum,  ne  =  1017  cm'3,  nK  =  6xl(T19  cm2/W,  I L  =  1014  W/cm2,  £  =  S,  and 
R0  =  260  fj. m ,  we  find  that  the  plasma  wavelength  Xp  =  2nd  cop  ~  106// m  and 
rg  *90. 

Taking  the  Fourier  transform  of  Eqs.(21a,  b)  together  with  Eq.(6)  in  the  variables 

oo 

z  and  r ,  with  the  convention  Q  =  (2tt)-1  Jq exp(ikz  -  ia>T)dkdco ,  gives 


In  obtaining  Eqs.(23)  we  have  set  n0  equal  to  unity.  In  the  following  section  Eqs.(23) 

and  (24)  are  solved  analytically  for  an  idealized  plasma  filament. 

c)  Analysis  of  EMP  in  Air 

We  consider  steady  state  and  non-steady  state  (transient)  propagation  of  the  laser 
pulse  for  an  idealized  case  that  can  be  solved  analytically.  The  laser  pulse  envelope  is 
taken  to  have  the  form 

|£i(0,z,r)|  =  El(z) sin(/Tr /tl)  ,  (25) 

for  0  <  t  <  rL  and  EL  =  0  otherwise.  The  front  of  the  laser  pulse  is  at  r  =  0  and  the 
back  is  at  r  =  tl  .  We  analyze  Eqs.(23)  and  (24)  for  the  case  of  an  idealized  plasma 
filament  with  radial  extent  rp  in  which  the  plasma  density  is  constant  for  r  <  rp  and  zero 
for  r  >  rp.  It  is  also  assumed  that  the  recombination  process  is  negligible  and  that  the 
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ionization  process  is  sufficiently  rapid  that  the  plasma  density  does  not  vary 
longitudinally  within  and  behind  the  laser  pulse,  i.e.,  cop  is  independent  of  r  .  Under 

these  assumptions,  Eqs.(23a)  and  (23b)  can  be  combined  to  give  equations  for  the  axial 
field  inside  (r  <rp)  and  outside  (r  >  rp)  the  plasma  filament  respectively, 


K2p{co,k) 


4  n 

co2{\  +  ive/(o)-co2po 


k]  (<y,  k)  Sz  (r,  co,  k)  -  (a  /  vg  +  k) -  (r  Sr  ( r,k,a )) 


,  (26a) 


\_d_(  A' 

rdry  dr j 


(26b) 


where 


2  r  n  co  2  cole  ,2 

K0(<B,k)  =  2  2  +  -  +  k2. 

VE 


(27a) 


K2p{(o,k)  =  K20{a>,k)  +  — 

C  (1  +  IV  e!  CO) 

and  oo po  denotes  the  constant  plasma  frequency  within  the  channel.  The  radial 


(27b) 


component  of  the  source  term,  i.e.,  Sr ,  can  be  neglected  provided  copo  rp/c»\.  This  is 

assumed  to  be  the  case  in  writing  Eqs.(26)  and  (27).  However,  in  the  numerical 
simulations  the  radial  source  term  is  not  necessarily  neglected. 

The  boundary  conditions  across  the  plasma-air  interface  at  r  =  rp ,  are  the 

continuity  of  the  axial  electric  field,  Ez{r=r~ )  =  Ez(r  =rp) ,  and  the  azimuthal 

magnetic  field  Be  {r = rp  )  =  Be(r=rp).  Applying  the  boundary  conditions,  the  axial 
electric  fields  inside  and  outside  the  plasma  filament  are  respectively, 


Ez(r,co,k )  = 


4?r  Sx{(o,k) 

<y2(l  +  ive/co)  - 


KorpKx(KorP) 

D(co,  k) 


to(*Pr) 


Ez{r,co,k) 


co2{\+ivJco)  D(o),k)  {kp)  M°n 


(28b) 
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where  I0  and  K0  are  the  modified  Bessel  functions  of  order  zero  of  the  first  and  second 
kind,  respectively, 

D(a),k)  =  h(co,k)KprpK0(K0rp)I,{Kprp)  +  K0rp  Kfx0rp) 1 0(tc prp) ,  (29a) 


and 


h{co,k )  = 


1  - 


co 


po 


co  (1  +  ivj co) 


K20{co,k) 


x(co,k) 


(29b) 


Making  use  of  Eqs.(28a)  and  (28b)  the  ratio  of  the  field  at  the  boundary  to  the  field  at  the 
center  of  the  filament  is 


EJf  =  rp)  _  hj^kyKpTp  Ii(jcprp)Ko(K0rp) 
Ez{r  =  0)  D{co,k)  -  K0rpKfic0rp) 


In  the  plasma  filament  there  is  a  radial  component  of  plasma  current  density  that  induces 
a  surface  charge  at  the  plasma-air  boundary.  The  resulting  charge  distribution  is  such 
that  the  EMP  fields  outside  the  filament  are  extremely  small.  Neglecting  the  radial 
source  term  Sr  in  Eq.(24),  continuity  of  the  azimuthal  magnetic  field  across  the 
boundary  implies  that 


E  (r=r+)  co 2 

r  v  p  /  j  ^  po 

Er  (r=rp)  co2(l  +  ivjco) 


(31) 


In  the  weakly  collisional  regime  (ve«  co)  and  for  co  ~  cop  it  follows  that  the  field 

outside  is  nearly  completely  shielded.  In  the  highly  collisional  regime,  however,  the 
EMP  field  is  not  shielded  outside  the  filament. 

In  subsections  ii)  and  iii)  we  analyze  and  discuss  limiting  cases  of  these  equations 
for  comparison  with  the  numerical  results. 
i)  Estimate  of  EMP  Power  Conversion  Efficiency 

A  rough  estimate  for  the  efficiency  of  converting  laser  power  to  EMP  power  can 
be  obtained  by  using  rather  general  arguments  concerning  the  EMP  field  amplitudes. 

This  provides  a  consistency  check  with  the  calculation  of  efficiency  obtained  by  the  far 
more  precise  analysis  which  follows.  In  addition,  this  estimate  shows  that  the  low 
efficiencies  obtained  by  more  precise  analysis  are  to  be  expected. 
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The  maximum  amplitude  of  the  longitudinal  electric  field  (wakefield)  induced  by 
the  driving  laser  pulse  propagating  in  the  plasma  filament  in  the  absence  of  collisions  is 
given  by  [  1 9] 


2  2  2 
_  n  me  aT 

E  « - — 

z,  max 


2  q  crL 


(32) 


where  aL  -  q\EL\/  mcco0  is  the  laser  strength  parameter.  The  frequency  of  the  wakefield 
is  in  the  THz  range  and  determined  by  the  laser  pulse  duration  rL .  In  the  absence  of 
collisions,  the  maximum  field  amplitude  occurs  when  the  laser  pulse  duration  is 
approximately  equal  to  the  plasma  frequency  cop  rL  ~  1 ,  [29].  The  other  components  of 

the  EMP  fields  induced  by  the  driving  laser  pulse  are  of  the  same  order  or  smaller  than 
the  above  maximum  axial  field.  The  power  in  the  EMP  is  estimated  to  be 

p EMP  ~  c ^2, max  2/8^  where  I  is  the  area  through  which  the  EMP  energy  flows.  The 
area  of  the  cylinder  enclosing  the  EMP  fields  is  estimated  to  be 
E  ~  7t{2NprpcTL  +  rp),  where  Np  is  the  number  of  periods  (i.e.,  plasma 

wavelengths)  of  the  wakefield  behind  the  laser  pulse.  The  parameter  Np  is  a  measure  of 
the  axial  extent  of  the  emitting  region  of  the  EMP.  The  driving  laser  pulse  power  can  be 
written  in  terms  of  the  laser  strength  parameter,  PL  =(c/\6)(mcco0  R/q)2  af .  Using 

Eq.(32)  and  the  expression  for  the  laser  power  we  find  that  an  estimate  for  the  EMP 
power  efficiency  is 


7  Power 


p 

1  EMP 

n1 

aL  A 0 

2 

r2NprpCTL  +  rp 

Pl 

~  T 

V  CtL  J 

R 2  J 

(33) 


where  X0  is  the  laser  wavelength.  As  an  illustration,  consider  a  laser  pulse  intensity  of 

IL  —  7  xlO13  W / cm2  and  wavelength  Xa  —  0.8// m,  the  laser  strength  parameter  for  these 

values  is  aL  =  5.7x1  O'3 .  For  a  laser  pulse  length  of  tL  =  250 / sec ,  the  power 
efficiency  is 


Vpo 


2x10 


-8 


2NprpCTL  +  rt 
R 2 


2A 


(34) 
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Taking  the  laser  spot  to  be  R  -  200pm ,  plasma  filament  spot  size  to  be  rp  =100 pm, 
number  of  plasma  wavelengths  behind  the  pulse  to  be  Np  ~  1 0 ,  the  efficiency  is 

extremely  small,  rjPower  ~  10~7.  This  estimate  is  typical  of  the  values  obtained  by 

numerically  solving  the  full  set  of  EMP  equations. 
it)  EMP  in  the  Steady  State  Regime 

In  the  steady  state  propagation  regime  it  is  assumed  that  both  the  plasma  density 
and  the  laser  pulse  are  independent  of  the  propagation  distance,  i.e., 

2  A  2 

d copo  fdz  =  8  El  Idz  -  0.  Even  in  a  homogeneous  medium  there  is  a  minimum 

propagation  distance  needed  to  reach  a  steady  state.  The  steady  state  EMP  propagates 
with  velocity  equal  to  the  phase  velocity  of  the  laser  pulse,  vg  .  The  transient  EMP  on  the 

other  hand  propagates  with  velocity  ~  d  n0.  The  propagation  distance  needed  to  reach  a 

steady  state  is  Lsteady_state  ~  rL  /|l  /  vg  -  n0  /c|  ~  2y1g  cx  L .  For  the  purpose  of  illustration  it 

will  be  assumed  that  the  collision  frequency  can  be  neglected. 

The  axial  EMP  electric  field  is  vanishingly  small  outside  of  the  plasma  filament, 
while  inside  the  filament  it  is  given  by 


~  4  7tSz(co) 

Ex  (r ,  co)  =  —  ^  - 


co  -co 


po 


1 


Io(Kpr) 

1 0  (K  p  rp)  J 


(35) 


The  EMP  magnetic  field  inside  and  outside  of  the  plasma  filament  is  given  by 

4 niSz{co)  IMpr) 


Be(r,co)  = 


®  I( AT„r  )’ 

p  O  V  p  P  J 


(36a) 


and 


CO  CKP  Ia (kp  rp )  Kx (fca  rp ) 

respectively.  For  typical  parameters  we  find  that  Kprp  ~  co  rp ! c  »  1  and 


(36b) 


KorP  ~  03 rP  l(cYg)  «  1  ar>d  the  EMP  magnetic  field  inside  and  outside  the  plasma 
filament  reduces  to 
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Be{r,co)  = 


4  ni  Sz{co) 


co 


2  ncr. 


V'2 


V  03 p°  J 


exp( -coporplc)Ix{coporlc),  (37a) 


Be{r,co)  = 


AniSz{co)  rp 


co. 


-Ki  {cor  !\  gyg) 


V  V 

po  g  ‘  g 

In  the  steady  state  and  in  the  absence  of  collisions  the  Fourier  transform  of  the  source 
term  in  Eq.(8)  is 


(37b) 


S2(a)  =  ~ 


iq  co  co 


2 me  c  co 


PO_  T 

2  lL 


{<o). 


(38) 


where  IL{ r)  =  c\El\2  / 8^  is  the  laser  intensity. 

To  obtain  the  temporal  variation  of  the  fields  it  is  necessary  to  invert  the 
corresponding  Fourier  transforms.  To  derive  an  analytic  expression  for  Be  (r,  r)  we 

assume  an  intensity  profile  IL  (r)  =  IL0  for  0  <t  <tl  and  zero  otherwise.  Substituting 
Eq.(38)  into  Eq.(36)  and  taking  the  inverse  Fourier  transform  results  in 

Be{r,z)  =  B0—H(r,t),  (39) 

r 

where 


H{r,r)  = 


((^/vgr,)2  +  r2)'2  {(r/ygrg)2  +  (r  -  rL)2) 


il/2 


(40a) 


and 


B 


o 


2nq  copo 

2  2  ^ 
me  6)0 


L0  * 


(40b) 


is  the  peak  magnetic  field  at  r  =  rp  and  r  =  rL.  Note  that  the  magnetic  field  in  Eq.(39) 

is  proportional  to  the  square  root  of  the  plasma  density.  This  dependence  is  due  to  the 
near-cancellation  of  the  axial  components  of  the  convective  and  inductive  current 
densities.  Evaluating  the  magnetic  field  at  r=rp  in  Eq.(37a)  and  noting  that 


0)Po  cp/c  »  1 ,  we  obtain  Be  (rp )  =  4  ni  Sz  {co)  l{copo  co) .  The  inverse  transform  of  this 
expression  leads  to  the  magnetic  field  given  in  Eq.(40b) 
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Figure  2  plots  the  distribution  of  the  axial  Poynting  flux,  i.e.,  radiance, 

I z  ( r ,  t)  =  cB2e  ( r ,  r)  /  An  in  the  ( r,  r )  plane  for  r>rp  for  laser  and  plasma  parameters 


rjo  lc  =  25,  tlcd  !2n  =  0.8,  cojco  =125,  *  =  8,  v  *140, 


X0  =  0.8  pm  (2p  =100  /um),  and  peak  laser  intensity  IL0  =  7xl013  W/cm2 .  For  these 
parameters  the  magnetic  field  at  the  surface  of  the  filament  is  B0  «  9G.  The  axial  flux  is 
localized  in  r  to  the  vicinity  of  the  laser  pulse  and  reaches  a  maximum  of  ~  1 8  kW/cm2 
at  mid-pulse  (r  =  rL  12 ).  In  the  radial  direction,  the  axial  flux  extends  many  plasma 
wavelength  beyond  the  plasma  column. 

The  EMP  power  propagating  in  the  axial  direction  resides  outside  of  the  plasma 
filament  and  is  given  by 


Pemp(t)  =  -^\Bl{r,T)2nrdr  =  B 2  r2p  A(r, xL ) , 
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(41) 


where  A =  ln[(l  +  f2)(l  +  Af2  ]  -  2sgn[fAf  ]  ln[ , 


2f 2  Af 2 +f2 +Af 2 +2  f  Af 

(f 

2  At 

2  a 2  .  A2  ,y/2 

+  T  +A  T  +lj 

f2  +Af 2  +2 

T 

<1 

’  f  =  vgygr/rp. 


tL  =  vg  ygTL  /  rp  and  Af  =  f  -  tl  .  Figure  3  plots  PEMP  (r)  normalized  to  the  peak 
laser  power  PL  for  the  same  parameters  as  in  Fig.  2.  The  power  ratio  is  typically  very 
small,  with  a  peak  value  of  ~4  x  10“!0 .  The  energy  in  the  EMP  pulse  is 


—  jdTPEMP(t)  —  rp  J<frA(r,ri),  (42) 


The  energy  efficiency,  defined  as  the  ratio  of  the  EMP  energy  to  the  laser  energy,  is  given 
by 


V  Ei 


w. 
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o  L 
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o  \ 


r_P^PQ 

Ro°>oJ 


(43) 


where  P0  =  m2c5 1 q2  =8.75GW, 


'¥(fL)  =  T?  f  dr  (44) 
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and  the  laser  pulse  energy  is  WL  =  (c/1 6)  EL  R20  rL  .  The  function  x¥(tl )  is  plotted  in 

Fig.  4  which  shows  that  it  is  a  slowly  increasing  function  of  the  variable  rL .  The  EMP 
conversion  efficiency,  Eq.(43),  is  plotted  in  Fig.  5.  For  typical  laser  parameters  the 
efficiency  is  extremely  small,  rj^  ~  4  x  10"'° . 

in)  EMP  in  the  Transient  Regime 

In  this  subsection  EMP  generation  for  non-steady  state,  i.e.,  transient,  pulse 
propagation  is  discussed.  To  understand  radiative  EMP  we  employ  a  heuristic  model  in 
which  the  non-steady  state  regime  is  represented  by  a  medium  varying  periodically  in  the 
pulse  propagation  direction.  This  can  represent  for  example  the  periodic  focusing  and 
defocusing  of  the  laser  pulse.  Since  all  the  field  components  can  be  expressed  in  terms  of 
the  axial  electric  field  we  spatially  modulate  this  axial  field  by  the  factor 
(l  +  £-cos(Kz)),  where  K  >  0  is  the  wavenumber  associated  with  the  non-steady  state 
component  of  the  ponderomotive  force  and  e  is  a  measure  of  the  modulation.  In  the 
spatially  modulated  medium  the  ponderomotive  force  has  superluminal  Fourier 
components  for  a  range  of  K.  For  these  values  of  K  the  EMP  is  radiative,  resulting  in  net 
power  flow  in  the  radial  direction.  This  can  be  understood  by  noting  that  the  Fourier 
transform  of  the  axial  electric  field  outside  of  the  filament  is  given  by  [see  Eq.(23a)], 

f  1  d  f  0  2/  V 

~JZ  "Tv  ~  K^K)  E^K">  =  °>  (45) 

\rar\  or)  ) 

where 

k2±(o),K)  =  k2(co,±K)  =  ±  2— K  +  K2  =  j  —  ±  K  1  -  (46) 

ygrg  vg  (v,  J  c2 

Following  Fourier  inversion  the  solutions  of  Eq.(45)  are  expressible  in  the  form 
K0  (k±  r )  exp(  ±  /  K  z )  exp(-/ cor ) .  These  solutions  represent  the  radial  variation  of  the 
forward  propagating,  exp(z  K  z) ,  and  the  backward  propagating,  exp(-z'Kz) , 
contributions  to  the  EMP,  respectively,  where  the  transverse  wavenumber  is  k±  . 

For  the  purpose  of  discussing  the  nature  of  the  solutions,  i.e.,  radiative  and  non- 
radiative  EMP,  we  consider  the  limit  of  zero  collision  frequency,  ve  =  0 .  For  k\>  0  the 
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solution  of  Eq.  (45)  represents  exponentially  decaying  non-radiative  solutions,  i.e., 

*0(l  K±  I  r)  ~  r~U 2  exP(~  |  k±  |  r) .  Note  that  the  exponentially  increasing  solution  is  non¬ 
physical.  Since  co  and  K  have  been  chosen  to  be  positive,  k+  is  always  positive  and 
represents  exponentially  decaying  solutions  that  are  non-radiative.  However,  k2  can  be 
less  than  zero  for  values  of  K  in  the  range 


2  co 


>  K  > 


co 


2  c  y 


2  • 


(47) 


In  this  case  k_  =  ±/|  k_  | and  the  solution  K0(jc_r)  for  large  r  represents  radiative 
solutions,  i.e.,  K0(K_r)  ~  r~U2  exp(  +  /|  k_  |r).  Here,  the  upper  sign  represents  a  radially 

incoming  wave  and  is  not  physical.  The  lower  sign,  however,  represents  a  radially 
outgoing  EMP  field  and  corresponds  to  radiative  fields.  The  source  of  the  radiative  field 
is  the  ponderomotive  force  containing  Fourier  components  with  phase  velocities  that  are 
superluminal.  In  the  laser  pulse  frame,  the  Cherenkov-like  emission  angle  of  the  EMP 
radiation  with  respect  to  the  laser  pulse  propagation  direction  is 

©c  =tan-1(|/r_|/A').  (48) 

The  phase  velocity  associated  with  the  Fourier  components  of  the  ponderomotive  force  is 
given  by 


co  v 


= 


g 


CO 


v«  K 


(49) 


and  the  EMP  is  radiative  when  vph  >  v g  .  The  radiative  and  non-radiative  branches  of 

the  EMP  can  be  illustrated  with  a  dispersion  diagram  shown  in  Fig.6.  The  frequency  and 
wavenumber  of  the  EMP  are  given  by  the  intersection  of  the  electromagnetic  mode  and 
the  vph  >  c  dispersion  line  shown  in  Fig.  6. 


III.  EMP  Generated  on  the  Surface  of  a  Dielectric 

In  this  section  we  analyze  the  generation  of  EMP  from  dielectric  surfaces.  The 
incident  laser  pulse  ionizes  the  dielectric  and  produces  an  oscillating  plasma  current  sheet 
which  is  a  source  of  EMP.  In  addition,  we  assume  that  the  dielectric  may  have  a  non¬ 
vanishing  second  order  susceptibility,  \  this  is  the  case  in  crystalline  dielectrics.  A 
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second  order  susceptibility  can  result  in  optical  rectification  and  be  a  source  of  EMP 
radiation  [43].  Other  methods  for  generating  THz  radiation  in  solids  are  discussed  in,  for 
example,  [22-25,  44], 

In  previous  sections  we  found  the  efficiency  of  laser  energy  conversion  into  EMP 
energy  in  air  to  be  extremely  small.  In  this  section  we  consider  the  EMP  generated  when 
a  laser  pulse  is  incident  on  a  dielectric  target.  A  difference  between  a  dielectric  and  the 
atmosphere  is  that  the  electron-neutral  collision  frequency  is  nearly-four-orders— of- 
magnitude  larger  in  dielectrics.  Another  difference  between  a  dielectric  and  the 
atmosphere  is  related  to  the  ionization  process.  Typically,  to  elevate  (ionize)  an  electron 
from  the  valence  band  into  the  conduction  band  in  a  dielectric  requires  an  amount  of 
energy  on  the  order  of  a  few  electron  volts.  Thus  for  a  laser  wavelength  of 
K  =0-8  fim  ( hco0  =1.54eV)  a  few  photons  suffice  to  induce  photoionization.  That  is, 


in  Eq.(5),  l  «  2  and  hence  the  plasma  filament  is  no  longer  limited  to  a  narrow  region 
near  the  center  of  the  laser  pulse,  i.e.,  the  filling  factor  is  nearly  unity.  Finally,  due  to  the 
short  penetration  depth  into  a  dielectric  [46],  self- focusing  and  diffraction  of  the  laser 
beam  are  negligible  and  consequently  a  ID  analysis  of  the  EMP  generation  suffices. 
Figure  7  is  a  schematic  illustrating  the  electromagnetic  fields  involved  in  our  analysis  of 
EMP  generation  when  an  ultrashort  laser  pulse  is  incident  on  a  dielectric. 

The  wave  equation  for  the  EMP  field  inside  the  dielectric  is 


VxVxE  + 


n l  d2E 
c 2  dt 2 


4  n  dJ 
c2  dt 


4n  d2  P2 
c2  dt 2  ’ 


(50) 


where  na  is  the  linear  index  of  the  dielectric,  the  plasma  current  density,  J ,  is  given  by 

I  ~  1 2 


Eq.(6),  P2  =  Xl  \El 


ex  /  2  is  the  second  order  polarization  field  polarized  in  the  x 


direction,  x2  is  the  second  order  susceptibility  and  EL(z,t )  is  the  amplitude  of  the  laser 
pulse  inside  the  dielectric  [41,43],  In  writing  the  second  order  polarization  field,  2nd 
harmonic  terms  have  been  neglected.  A  typical  value  for  the  second  order  susceptibility 
in  crystalline  dielectrics  is  x2  ~  10“7  esu.  On  the  right  hand  side  of  Eq.(50),  the  ratio  of 

the  second  order  polarization  field  source  term  to  the  radial  and  axial  components  of  the 
plasma  current  density  source  term,  are  respectively 
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|(dP2/dO/J-er| 

„  mcco]  x2(2  veR0^ 

(51a) 

4  (o\  tl  \  c  y 

|(3P2  /  dt)l  J  -e2|  - 

1-  mCCOl  ^2 

(51b) 

^  rL 

In  obtaining  Eqs.(51)  it  was  assumed  that  the  process  is  in  the  collisionally  dominant 
regime,  ve  tL  »  I .  For  coQl  cop  ~  10,  rL  ~  lOOfsec,  ve  ~  1015  sec-1  ,  and  spot  size  of 

Ra  ~  200pm  the  ratios  in  Eqs. (51a, b)  are  ~105  and  102  respectively.  This  indicates 
that  the  optical  rectification  process,  compared  to  the  plasma  generation  process,  can  be  a 
far  more  efficient  mechanism  for  generating  EMP  in  crystalline  materials. 

To  obtain  both  the  current  density  and  second  order  polarization  field  on  the  right 
hand  side  of  Eq.(50),  the  laser  pulse  amplitude  EL  is  required.  To  obtain  the  laser  pulse 
amplitude  we  write  the  incident,  reflected  and  transmitted  laser  pulse  fields  as  follows 
E,„(z>0  =  (1/2) Ein (t-z/ c) exp[/ (k0 z  -  a>0t)]tx  +c.c.,  (52a) 

ER(z,t)  =  (l/2)ER(t +  z  /  c)exp[-i(k0z  +  coot)]ex  +  c.c.,  (52b) 

EL(z,t)  =  (l/2)EL(z,t)exp(-ico0t)ex  +  c.c.,  (52c) 

where  k0  -  co0l  c  in  the  vacuum  region.  The  laser  field  inside  the  dielectric  satisfies 

a2E,  1  82EL  An  dJ,  An  d2  P 

=  +  <53) 

where  J L(z,t)  and  P(z,/)  are  respectively,  the  current  density  and  total  polarization 

fields.  Within  the  dielectric  the  current  density  is  given  by  Eq.(2) 

where 

JL(z,t)  =  (l/2)JL(z,t)exp(-ia)0t)ex  +  c.c.,  (54) 

is  the  laser  driven  plasma  current  density.  The  rate  of  change  of  plasma  density  is  given 
by  Eq.(4a).  The  polarization  field  consists  of  a  linear  and  nonlinear  contribution, 

P (z,0  =  (1/2 )(PL(z,t)  +  PNL(z,t))txp(-icoo  t)ex  +  c.c.,  (55) 

/S  A  A  |  A  2  a 

where  PL(z,t)  =  %L  EL(z,t) ,  PNL  =  a\EL(z,t)  EL(z,t),  Xl  is  the  linear  susceptibility 

and  a  is  related  to  the  nonlinear  Kerr  refractive  index.  Note  that  the  second  order 
polarization  field  in  Eq.(55)  is  omitted  since  it  does  not  contribute  to  a  driving  term  at  the 
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laser  frequency.  If  the  laser  pulse  duration  is  long  compared  to  the  laser  period,  i.e., 
co 0  »  \dl  dt\ ,  Eqs.(53)  and  (2)  can  be  written  as 

a 
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JL(z,t)  = 


-®2p(z,t)EL(z,t), 


An{co0  +  ive) 

where  nQ  =  (1  +  An  %L)h2  is  the  linear  refractive  index.  To  lowest  order  in  81  dt , 
Eqs.(56a)  and  (56b)  combine  to  give 
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V  lMP  ) 


and  con  =Anq  n,l0/m  ~  1.8xl016sec  1  for  «  ~  1023  cm  3.  Note  that 


Ana 


El  ( z,t)\  =2  n0nK  («„  /  nn0 )  I L  (z,  t)  where  is  the  nonlinear  Kerr  refractive 


index,  «(J  -  nno  -  ne  and  IL{z,t)  =  cn0  EL(z,t)\  !%n  is  the  laser  pulse  intensity.  Fora 

dielectric  the  number  of  photons  needed  for  ionization  is  typically  small,  i.e., 
f  =  1,  2,  or  3,  for  Ao  =  0.8pm  .  Equation  (57)  together  with  Eqs.(58)  describe  the 


evolution  of  a  laser  pulse  propagating  in  an  ionizing  dielectric.  In  dielectrics  collisional 
and  ionization  absorption  of  the  laser  pulse  takes  place  over  a  short  distance,  i.e.,  the 
penetration  depth  is  a  few  microns.  For  distances  short  compared  to  the  laser  pulse 
length,  i.e.,  z  «  crL  / n0,  transit  time  effects  can  be  neglected.  In  this  limit  8/dt  can  be 
neglected  and  Eq.(57)  reduces  to 
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+  K2(z,t)  E,(z,t)  =  0. 


(59) 


The  WKB  solution  of  Eq.(59)  is 

(K( 0  t)  Y/2  f  r 

^(z,0  =  £,(0,0  -^4  exp  /jVz'£(z',0 
\K{z,t)j  y  0 

where  solutions  representing  waves  traveling  towards  the  left  are  not  considered.  The 
value  of  the  field  at  the  boundary,  i.e.,  EL  (0,0 ,  is  obtained  from  the  boundary 
conditions.  Applying  the  boundary  conditions  at  z  =  0 ,  yields 

Ein (t)  +  ER(t)  =  £,(0,0  and 

Ein (/)  -  ER(t)  =  EL(0,t)(K(0,t)+(i/2)(d\n(K)/dz)2=0)/ k0  which  combine  to  give 


£,(  o,r)  = 


2  KEJt) 

,  ^/A  ,  i(d\n(K) 

K  +  K(0,t)  +  -  — ^ 
21  oz 


Within  the  dielectric  the  laser  pulse  amplitude  is  given  by  Eq.(60)  together  with  Eq.(62). 
The  laser  pulse  amplitude  determines  the  current  density  given  by  Eq.(2)  and  the  second 
order  polarization  field  which  drives  the  EMP  fields  in  Eq.(50).  These  equations  are 
solved  numerically  in  Section  IV. 


IV.  Numerical  Simulations  of  EMP  Generation 

The  numerical  simulations  presented  in  this  section  are  based  on  solving  Eq.  (20) 
for  the  EMP  field,  with  the  plasma  current  and  source  given  by  Eqs.  (6),  (8)  and  (10). 

The  laser  envelope  evolves  according  to  Eqs.  (14)  and  the  ionization  generated  plasma 
density  is  given  by  Eq.(l  5c).  The  characteristic  ionization  and  recombination  parameters 
used  in  the  simulations  are  discussed  following  Eq.  (5).  The  system  of  equations  is 
solved  in  axisymetric  cylindrical  geometry  and  evolved  using  a  time  and  space  centered 
finite  differencing  scheme  with  the  electric  field  and  plasma  density  calculated  on  the 
same  grid  and  the  magnetic  field  and  plasma  current  calculated  on  a  separate,  staggered 
grid.  The  equations  for  the  EMP  fields  are  solved  in  the  laboratory  frame  characterized 
by  the  variables  ( r,z,t ) .  The  propagation  equations  for  the  driving  laser  pulse  are  solved 


26 


in  the  group  velocity  frame  and  transformed  to  the  laboratory  frame  when  coupled  to  the 
EMP  equations. 

The  SDE  method  [42]  of  solving  for  the  laser  pulse  amplitude  provides  a 
description  of  the  laser  pulse  evolution  that  is  simple  to  implement  numerically. 
However,  it  is  based  on  the  premise  that  the  transverse  laser  profile  is  approximately 
Gaussian.  In  certain  regimes  this  may  be  violated  as,  for  example  in  Fig.  13  of  Ref  [36], 
where  a  hollowing  of  the  laser  pulse  takes  place.  In  such  cases,  the  solution  of  Eq.  (1 1) 
must  be  employed.  In  the  following  we  shall  assume  that  the  Gaussian  profile 
approximation  is  valid. 

The  simulations  are  initialized  at  t  =  0  with  the  laser  pulse  propagating  in 
vacuum  in  the  positive  z  direction.  The  leading  edge  of  the  laser  pulse  at  t  =  0  is  located 
at  axial  position  z  =  0 .  In  the  vacuum-air  transition  region  from  z  =  0  to  z  =  Ln ,  the 
neutral  density  increases  with  the  functional  form 

„  (2)=„  f  Z>-L • 

*“[[10(z/£„)!  -  15 (z/£„)‘  +  6(z/I„)s],  z  S  L, 
and  has  the  constant  value  nllo  =  2.7  x  1019  cm-3  for  z  >  Ln. 

a)  EMP  Generation  in  Air 

The  following  four  examples  illustrate  the  generation  of  EMP  from  the  interaction 
of  an  ultrashort  laser  pulse  with  air.  The  first  two  examples  validate  the  numerical 
simulation  with  the  analyses  of  Sec.  III.  The  third  example  is  a  simulation  of  a  recent 
experiment  which  observed  EMP  emitted  from  a  laser  pulse  propagating  in  air.  The  last 
example  in  this  subsection  shows  EMP  generation  from  a  much  higher  intensity  laser 
pulse  which  fully  ionizes  the  air. 
i)  Steady  State  Laser  Pulse  Propagation 

This  first  example  benchmarks  the  numerical  simulation  with  the  steady  state 
analysis  of  Sec.  Ilia.  A  laser  pulse  propagates  through  a  preformed  plasma  filament 
characterized  by  a  constant  density  ne  =10,7cm'3  (<x>po  =1.9xl013  sec'1)  in  the  region 

r  <rp  and  ne  =  0  for  r  >  rp ,  where  the  filament  radius  is  rp  =  0.4  mm  (~  25  c/copo) . 
The  laser  spot  size  is  assumed  to  be  much  larger  than  rp  so  that  the  transverse  variation 
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of  the  laser  envelope  can  be  neglected.  To  approximate  the  square  laser  pulse  profile 
assumed  in  the  analysis,  the  laser  envelope  in  the  simulation  is  taken  to  have  a 
longitudinal  profile  given  by  El(t )  =  EL0  sin2  [  n(r/rL) 16 12]  for  0  <r  <rL  and 
El  (r)  =  0 ,  otherwise.  The  smooth  variation  of  the  laser  envelope  at  the  leading  and 
trailing  edges  reduces  numerical  dispersion.  The  pulse  duration, 
rL  =  270  fsec  ~  0.8  x  2 n  I copo ,  is  chosen  to  optimally  excite  a  plasma  wave  under 

collisionless  conditions.  The  peak  laser  intensity  is  IL0  =  7xlOn  W/cm2  and  the  laser 
wavelength  is  X0  =  0.8  pm .  To  be  consistent  with  the  steady  state  analysis,  the  laser 
envelope  is  independent  of  the  propagation  distance  t)  and  the  radial  EMP  source  term 
given  by  Eq.  (10)  is  neglected  in  this  simulation  example. 

The  length  of  the  vacuum-air  transition  region  for  this  simulation  example  is 
Ln  =  20  ct  a>po  -  0.3  mm .  A  transient  EMP  is  generated  in  this  transition  region  which 

must  be  allowed  to  propagate  away  from  the  laser  pulse  before  the  steady  state  is 
achieved.  To  achieve  steady  state,  the  laser  pulse  must  propagate  a  distance 
z  »  2 y\ciL  .  To  reduce  computer  run  time  and  expedite  reaching  the  steady  state,  a 

vacuum  ( n0  =  1 )  is  assumed  outside  of  the  plasma  filament  and  the  laser  pulse  is  given 
an  artificially  low  group  velocity,  yg  =13.  For  these  parameters,  the  minimum 
propagation  distance  required  to  reach  steady  state  is  ~  3  cm . 

Figure  8  plots  the  longitudinal  laser  intensity  profile  and  axially  directed  EMP 
field  after  a  propagation  distance  of  z  =  6.3  cm  .  Solid  curves  denote  simulation  results 
and  dashed  curve  denote  the  results  of  the  steady  state  analysis  given  by  Eq.  (35).  Figure 
8a  compares  the  laser  intensity  profile  used  in  the  simulation  with  the  theoretical  square 
profile  used  in  the  analysis.  Figure  8b  plots  the  normalized  axial  EMP  field 

Ez  =  \q\EJ  mca)0  on  axis  as  a  function  of  normalized  time  r  Q)p0  /  In  and  Fig.  8c 
shows  the  radial  profile  of  Ez  at  r  copo  /  2zr  =  1.4.  The  frequency  of  Ez  behind  the  laser 

pulse  is  equal  to  the  plasma  frequency.  There  is  good  agreement  between  the  simulation 
and  the  theoretical  predictions  for  both  the  axial  and  radial  profiles  of  Ez .  The  slight 
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discrepancy  in  the  EMP  field  amplitude  between  the  simulation  and  the  analysis  apparent 
in  Figs.  8b  and  8c  is  attributable  to  the  difference  in  pulse  shape  shown  in  Fig.  8a. 

Figure  9  plots  the  normalized  EMP  magnetic  field  Be  =  \q\  Belmca>0  in  the 
steady  state  corresponding  to  the  simulation  of  Fig.  8.  Figures  9a  and  9b  plot  contours  of 
Be  in  the  (r,  r)  plane  obtained  from  the  steady  state  analysis  and  the  simulation, 
respectively.  The  plots  for  the  steady  state  analysis  were  obtained  from  Eq.  (39).  Figure 
9c  shows  radial  profiles  of  Be  taken  at  mid-pulse  ( r  =  rL  /  2 )  from  the  simulation  (solid 
curve)  and  the  analysis  (dashed  curve).  There  is  good  agreement  between  theory  and 
simulation  for  the  amplitude  and  spatial  distribution  of  Be  in  the  steady  state  regime. 

The  radial  EMP  Poynting  flux  Ir  =  cEzBg  /  An  (solid  curve)  and  laser  intensity 
profile  (dashed  curve)  outside  the  plasma  filament  is  plotted  in  Fig.  10a  as  a  function  of 
normalized  time  t  copo  /  2tt  .  Note  the  adjacent  regions  of  positive  and  negative  radial 

flux  which,  when  integrated  over  time,  nearly  cancel.  This  result  is  consistent  with  the 
analytical  prediction  that  the  EMP  field  is  non-radiative  in  the  steady  state  and  there  is  no 
net  radial  energy  flux.  Figure  10b  shows  that  there  is,  however,  a  much  larger  net  axial 
Poynting  flux  (~  400  W/cm2 )  in  the  simulation  which  is  in  agreement  with  the  analysis. 
The  axial  EMP  flux  is  localized  longitudinally  roughly  within  the  laser  pulse. 
ii)  Transient  Laser  Pulse  Propagation 

In  the  next  example,  we  verify  the  condition  for  radiative  EMP  given  by  Eq.  (47). 
A  laser  pulse  with  longitudinal  profile  given  by  EL  (r)  =  EL0  sin[  n  r  /  2  xL  ]  for  0  <  r  <  rL 

and  El(t)  =  0,  otherwise,  pulse  duration  rL  =  270  fsec  ~  0.8x2;r/typo,  peak  laser 

intensity  IL0  =  7  x  1013  W/cm2 ,  and  wavelength  A0  =  0.8  pm  propagates  through  a 

modulated,  preformed  plasma  filament.  The  plasma  density  is  given  by 

ne{r,z)  =  neo[  1  -  ©(r-^)]^  +  es\n{Kz)},  where  neo  =  1017  cm"3,  the  modulation 

amplitude  is  s  =  0.3  and  0(x)  is  the  Heaviside  unit  step  function.  As  in  the  previous 

example,  there  is  a  transition  region  from  vacuum  to  air  with  a  scale  length 

Ln  =20  d  co po  =  0.3  mm .  The  linear  and  nonlinear  (Kerr)  indicies  of  air  are  chosen  to 

be  7i o  —  1  =  2.7x10  and  nK  =6x10  19  cm 2 /W ,  respectively.  For  these  parameters,  the 
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relativistic  factor  associated  with  the  laser  group  velocity  [see  Eq.  (22)]  has  the  value 
y  =  109.  The  laser  envelope  is  taken  to  be  independent  of  the  propagation  distance  7 

and  the  radial  EMP  source  term  given  by  Eq.  (10)  is  neglected  to  be  consistent  with  the 
analysis.  Note  that  for  an  unmodulated  plasma  filament,  the  distance  required  to  reach 
steady  state  for  these  parameters  is  ~  2  m . 

Figures  1 1  and  12  compare  the  spatial  EMP  distribution  generated  from  two 
plasma  filaments  modulated  with  different  wavelengths.  In  Fig.  1 1,  the  modulation 
wavenumber,  K  =  copo  I  Ac ,  is  chosen  to  satisfy  the  condition  for  radiative  EMP  given  by 

Eq.  (47)  assuming  that  the  characteristic  frequency  co  ~  copo .  Figure  1  la  plots  the  spatial 

(g  =  z  -  vgt,  r )  distribution  of  the  radial  Poynting  flux,  Ir ,  after  a  propagation  distance 

of  ~2cm .  Radiative  EMP  is  observed  propagating  from  the  plasma  channel  behind  the 
laser  pulse.  The  radiation  is  confined  to  a  cone  characterized  by  the  Cherenkov  angle 
given  by  Eq.  (48).  The  outward  EMP  flux  observed  within  the  laser  pulse 
(-80  pm  <  £  <  0)  is  a  manifestation  of  the  transient  interaction  at  the  vacuum-air 
transition  region.  It  is  expected  that  there  would  be  no  net  outward  flux  in  the  region 
-  80  pm  <  g  <  0 ,  had  the  simulation  reached  the  steady  state  regime  after  2  m  of 

propagation.  In  Fig.  1  lb,  the  radial  profile  of  I r  behind  the  laser  pulse  at  axial  position 
=  -  220  pm  shows  an  outward  Poynting  flux  on  the  order  of  100  W/cm2  extending 
outside  the  plasma  filament. 

In  Fig.  12,  simulation  parameters  are  identical  to  those  of  Fig.  11,  except  that  the 
modulation  wavenumber,  K  =  Acopo  I  c ,  is  chosen  to  be  larger  than  the  maximum  value 

for  radiative  EMP,  assuming  co  ~  copo .  Consistent  with  the  analysis,  Fig.  12a  shows  an 

absence  of  radial  flux  behind  the  laser  pulse  outside  the  plasma  filament.  Behind  the 
laser  pulse,  the  EMP  is  confined  within  the  plasma  filament.  The  radial  profile  of  Ir  in 
Fig.  12b  shows  that  while  there  is  a  larger  radial  EMP  flux  within  the  plasma  filament 
relative  to  the  simulation  of  Fig.  1 1,  it  is  highly  evanescent  outside  of  the  plasma 
filament.  Again,  the  outward  flux  observed  in  the  region  -80  pm  <  <  0  is  due  to  the 

transient  interaction  at  the  vacuum-air  transition  region,  and  is  not  present  in  the  steady 
state  regime. 
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iii)  Comparison  with  Experiment 

A  recent  experiment  reported  in  Ref.  [17],  measured  an  EMP  from  a  femtosecond 
laser  pulse  propagating  in  air.  In  the  experiment,  a  laser  pulse  with  a  FWHM  duration  of 
120  fsec,  wavelength  Xo  -  0.8  pm ,  and  energy  ~30  mJ  was  passed  through  an  aperture 

of  diameter  1 .5  cm  and  focused  into  the  air  with  a  2  m  focal  length  lens.  As  the  laser 
pulse  focused  and  the  intensity  increased,  a  plasma  filament  was  formed  near  the  focal 
region  and  EMP  was  measured  along  the  length  of  the  plasma  filament.  The  following 
results  are  from  a  simulation  with  similar  parameters. 

In  our  simulation,  the  parameters  for  air  are  taken  to  be  n0  -  1  =  2.7  x  1 0  4 , 
nK=  3  x  10  '9  cm2/W,  Uion  =  12  eV,  £  =  8,  lMP  =  1014  W/cm2,  and 
nm>  =  2-7  x  1019  cm  3  .  The  laser  pulse  envelope  is  taken  to  be  Gaussian  in  the  transverse 
direction  with  a  spot  size  R0  =  7.5  mm .  The  longitudinal  laser  field  profile  is  given  by 
Eq.  (25)  with  rL  =  240  fsec ,  which  corresponds  to  a  FWHM  intensity  duration  of  120 
fsec.  The  laser  pulse  has  wavelength  X0  -  0.8  pm  ,  peak  power  PL  -  83  GW  and  energy 

~  10  mJ.  The  simulation  is  performed  in  two  parts  as  illustrated  in  Fig.  13.  The  first  part 
models  the  focusing  of  the  initially  low  intensity  laser  pulse  in  air  by  a  lens  with  a  focal 
length  of  2  m .  Here  the  simulation  covers  a  range  0  <  z  <  z0,  where  zo  =  1 .82  m . 

Over  this  range,  the  laser  intensity  is  sufficiently  low  that  no  appreciable  plasma  density 
is  formed  and  no  EMP  is  generated.  To  expedite  the  simulation,  only  the  propagation 
equations  for  the  driving  laser  pulse  envelope  are  solved  for  z  <  za.  The  laser  pulse 

envelope  at  z  =  z  is  used  to  initialize  the  full-scale  EMP  simulation  for  z  >  z  . 

o 

Figure  14a  plots  the  peak  plasma  density  (solid  curve)  and  laser  intensity  (dashed 
curve)  over  the  range  0  <  Az  <  16  cm ,  where  Az  =  z  -  z0  is  the  distance  from  the  start 
of  the  full-scale  EMP  simulation.  For  Az  <  9  cm,  the  laser  intensity  increases  due  to 
both  geometric  and  nonlinear  focusing.  For  Az  >  9  cm ,  a  plasma  filament  is  formed  on 
axis  which  defocuses  the  trailing  edge  of  the  laser  pulse,  thereby  resulting  in  a  shorter 
duration  laser  pulse.  The  peak  plasma  density  achieved  in  the  simulation  is 
~  2  x  1 0  cm  3 ,  which  occurs  at  Az  ~  1 0  cm .  For  Az  >10  cm  the  plasma  density 
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steadily  decreases  due  to  the  shortening  of  the  laser  pulse.  Figure  14b  plots  the 
corresponding  peak  values  of  the  axial  (solid  curve)  and  radial  (dashed  curve)  Poynting 
fluxes  as  a  function  of  A z .  The  maximum  EMP  flux  coincides  with  the  maximum 
plasma  density.  The  maximum  radial  flux  is  Ir  ~  35  kW/cm2  while  the  maximum  axial 
flux  is  Iz  ~  15  kW/cm2 .  Significant  shortening  of  the  driving  laser  pulse  prevented  the 
simulation  from  extending  beyond  Az  ~  16  cm  without  violating  the  assumption  that 
group  velocity  dispersion  is  negligible. 

Figure  1 5  plots  the  spatial  distribution  of  the  laser  intensity  and  the  various  EMP 
field  components  at  Az  ~  1 0  cm ,  when  the  plasma  density  is  maximum.  The  plasma 

filament  radius  in  these  plots  is  rp  ~  100  jjm .  The  hour-glass  shaped  laser  intensity 

contours  show  nonlinear  Kerr  self-focusing,  i.e.,  the  higher  intensity  middle  region  of  the 
pulse  has  a  smaller  spot  size  than  the  lower  intensity  leading  and  trailing  edges.  In 
addition,  plasma  defocusing  is  evident  from  the  slightly  larger  spot  size  in  the  trailing 
regions  of  the  pulse  where  the  plasma  density  is  higher.  The  axial  EMP  electric  field  is 
localized  within  the  plasma  filament  while  the  magnetic  field  and  radial  electric  field 
extend  beyond  the  transverse  dimension  of  the  filament. 

Figures  16a  and  16b  plot  the  spatial  distribution  of  the  radial  and  axial  Poynting 
flux  at  Az  ~  10  cm ,  respectively.  Inside  the  plasma  filament,  the  Poynting  fluxes  are 
peaked  within  the  laser  pulse  and  decay  in  amplitude  behind  the  pulse.  The  radial  flux  is 
mainly  confined  within  the  plasma  filament  while  the  axial  flux  extends  radially  beyond 
the  filament.  Although  not  apparent  from  the  color  scale  of  the  plot,  a  small  radial  flux  is 
present  outside  of  the  filament  which  can  be  experimentally  detected.  Figure  17a  plots 
the  on  axis  laser  intensity  (normalized  to  its  peak  value  at  z  =  0 )  versus  time  at 
Az  =  10  cm  (dashed  curve)  and  Az  =  12  cm  (solid  curve).  The  relatively  shorter  pulse 
duration  at  Az  =  12  cm  is  the  result  of  plasma  defocusing  of  the  trailing  edge  of  the  laser. 
Figure  17b  plots  the  corresponding  radial  Poynting  flux  versus  time  at  radial  position 
r  -  225  pm  .  The  Poynting  flux  is  peaked  in  a  region  behind  the  driving  laser  pulse  and 

has  a  duration  that  is  of  the  order  ~\l  ve.  Note  that  the  flux  at  Az  =  10  cm  is  directed 
radially  outward  while  at  Az  =  12  cm  ,  it  is  outward  just  behind  the  laser  pulse  and  then 
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reverses  direction.  Figure  17c  plots  the  frequency  distribution  of  the  two  curves  shown  in 
Fig.  17b.  The  frequency  distribution 


A  00 

s  -j==  jlr(r0,t)e“s‘  dt, 


(63) 


has  units  of  intensity  and  represents  the  intensity  contained  within  a  frequency  band  A co , 
which,  for  a  laboratory  experiment,  is  typically  the  resolution  of  the  detector  used  to 
measure  the  EMP .  In  this  case,  Aco I2n  ~  1  THz  is  the  smallest  frequency  interval  that 
can  be  resolved  by  the  simulation  grid  size.  The  EMP  spectrum  at  Az  =  10  cm  is  peaked 
near  zero  frequency  and  is  broadband  with  a  width  of  ~10  THz .  At  Az  =  10  cm ,  the 
maximum  electron  collision  frequency  ve  ~  6x  1013  sec  1 .  At  Az  =  12  cm ,  the  spectrum 


is  broader  and  a  local  maximum  at  a  frequency  of  ~  2  THz  is  apparent.  The  peak  plasma 
frequency  varies  from  ~  13  THz  to  ~  12  THz  over  the  2  cm  of  propagation  shown  in  the 
plots.  Note  that  there  is  no  distinct  feature  in  the  EMP  spectrum  at  the  plasma  frequency 
or  any  of  its  harmonics. 

Figure  1 8  plots  the  total  EMP  power  emitted  radially  through  a  cylindrical  surface 
of  radius  r0  =  580  pm  versus  time,  i.e., 


oo 

P(t)  =  2/T ro  | /,. (rQ , z, t)  dz  .  (64) 

-oo 

Note  that  time  t  =  0  in  the  plot  corresponds  to  the  instant  the  laser  pulse  reaches  the  axial 
position  z  =  z0 .  The  power  has  a  maximum  of  ~  0.8  W.  Most  of  the  power  is  emitted 

over  a  duration  ~  0. 1  nsec.  Integrating  the  curve  in  Fig.  1 8  over  time,  the  total  energy 
radiated  in  the  radial  direction  during  the  entire  simulation  is  found  to  be  ~  50  pJ. 
iv)  High  Energy  Laser  Pulse 

As  a  final  example  of  atmospheric  EMP  generation,  we  consider  the  propagation 
of  a  high-energy  (~  1  J),  high-intensity  ( >  1 014  W/cm2 )  ultrashort  laser  pulse  in  air.  In 
this  example,  the  initial  laser  pulse  is  Gaussian  in  the  transverse  direction  with  a  spot  size 
R0  =  1  mm  •  The  longitudinal  envelope  is  given  by 

|£i(0,z,r)|  =  EL(z)sin2(7TT/TL)  (65) 
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for  0  <  t  <  tl  and  zero  otherwise,  with  rL  =  240  fsec .  The  laser  pulse  is  initially 
collimated  and  has  wavelength  X0  =  0.8  pm  ,  peak  power  PL  =  5.3  TW,  intensity 
IL0  =3.5x1 014  W/cm2 ,  and  energy  ~  0.5  J.  Parameters  for  air  are  taken  to  be 
n0  -  1  =  2.7xl0'4 ,  nK  =  3  x  10'19  cm2/W,  U i0H  =  12  eV,  l  =  8,  and 
IMP  =  1014  W/cm2 .  The  air-vacuum  interface  has  scale  length  ~  70  pm  . 

Figure  19a  shows  the  normalized  peak  laser  intensity,  I L  / 1 L0  and  normalized 
peak  plasma  density  ne  /  nll0  as  a  function  of  propagation  distance.  At  z  =  0 ,  the  laser 

pulse  is  in  vacuum.  As  the  laser  pulse  enters  the  air,  the  air  is  fully  ionized  and  the  laser 
intensity  decreases  rapidly.  Note  that  for  the  initial  parameters  of  the  simulation,  the 
scale  length  for  laser  energy  depletion  due  to  ionization  is  ~  8  jum.  Over  a  propagation 
distance  of  ~  1  cm,  the  laser  pulse  is  seen  to  lose  more  than  ~  70%  of  its  energy.  As  the 
laser  intensity  decreases,  the  plasma  density  is  also  seen  to  decrease.  The  corresponding 
values  of  peak  radial  and  axial  Poynting  flux  are  plotted  in  Fig.  19b.  Initially,  the  peak 
radial  flux  (~  1  MW/cm2 )  is  much  larger  than  the  peak  axial  flux  (~  26  kW/cm2) . 

After  ~  1  cm  of  propagation,  the  radial  and  axial  flux  are  both  of  order  ~  1  kW/cm2 . 
Note  that  a  quasi-equilibrium  state  is  attained  after  ~0.3  cm  of  propagation  in  which  the 
laser  intensity  and  plasma  density  change  slowly  relative  to  the  initial  phase  of 
propagation  in  which  the  air  was  close  to  fully  ionized. 

Figure  20  plots  the  corresponding  spatial  distributions  of  laser  intensity,  and 
Poynting  flux  in  the  quasi-equilibrium  configuration  at  a  propagation  distance 
z  =  0.8  cm  .  The  radius  of  the  plasma  filament  is  these  plots  is  ~  0.5  mm .  Due  to  the 
high  electron  collision  frequency,  there  is  little  axial  flux  in  the  region  r  <  0.25  mm, 
while  the  radial  flux  is  seen  to  be  strongly  damped  behind  the  laser  pulse, 
b)  EMP  Generation  on  the  Surface  of  a  Dielectric 

In  this  subsection,  we  simulate  the  EMP  generation  from  the  interaction  of  a  laser 
pulse  with  a  dielectric  surface.  The  simulation  is  set  up  with  air  in  the  region  z  <  0 ,  and 
a  dielectric  in  the  region  z  >  0 ,  as  indicated  in  Fig.  7.  The  initial  laser  pulse  at  /  =  0  is 
situated  just  outside  the  dielectric  and  has  a  field  envelope  that  is  Gaussian  in  the 
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transverse  direction  with  a  spot  size  R0  =  0.1  mm.  The  longitudinal  laser  field  profile  is 
given  by  Eq.(25)  with  rL  =  200  fsec .  The  laser  pulse  has  wavelength 
X0  =  0.8  \xm(  ficoa  =  1.54  eV) ,  and  peak  intensity  IL  =  1013  W/cm2 .  The  dielectric  is 
characterized  by  nno  =  6  x  1022  cm'3,  n0  =3,  nK  =  1  x  10'15  cm2/W,  Uion  =3  eV, 

1  =  2,  and  IMP  =  5 x  1014  W/cm2 .  The  WKB  solution,  i.e.,  Eq.  (60),  is  evaluated 

numerically  to  obtain  evolution  of  the  laser  pulse  within  the  dielectric,  and  Eq.(50)  is 
solved  to  determine  the  EMP  generation. 

Figure  21  plots  contours  of  the  laser  intensity,  I L ,  and  EMP  intensity, 
l  emp  =C +  I?)1'2  in  the  (z,r)  plane  at  various  times  during  the  simulation.  Fig  21a 
shows  the  laser  intensity  contours  before  the  interaction  at  t  =  0.  Figure  21b  shows  the 
simulation  results  immediately  after  the  front  of  the  laser  pulse  enters  the  dielectric  and 
the  initial  EMP  is  generated.  The  EMP  extends  radially  on  the  order  of  the  laser  spot  size 
and  longitudinally  ~  50  pm  into  the  dielectric.  The  EMP  has  a  peak  intensity  of 

I emp  ~  30  MW/cm2  near  the  boundary.  At  this  time,  a  plasma  layer  of  width  ~5pm , 

radial  extent  ~  100  pm,  and  peak  density  ne  ~  3.5xl021  cm'3  ( cop  ~  3.3xl015  sec"1)  is 

formed  at  the  boundary.  The  laser  pulse  is  absorbed  within  several  microns  of  the 
boundary. 

Figure  21c  shows  the  EMP  as  the  trailing  end  of  the  laser  pulse  approaches  the 
boundary.  At  this  instant,  the  peak  EMP  intensity  is  I EMP  ~  1  MW/cm2  within  the 
dielectric.  Outside  the  dielectric,  an  EMP  traveling  in  the  negative  z  direction  is 
observed.  Figure  2 Id  shows  simulation  results  later  in  time  after  the  laser  pulse  has  been 
completely  absorbed  in  the  dielectric.  There  is  relatively  little  EMP  flux  within  the 
dielectric;  most  of  it  has  been  either  reflected  or  absorbed  by  the  plasma.  At  this  time, 
the  peak  plasma  density  is  ne  ~  9  x  1 021  cm  3 .  Within  the  dielectric,  a  pulse  of  intensity 
~  300  kW/cm2  traveling  in  the  positive  z  direction  is  observed. 
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V.  Summary 

In  this  paper  the  generation  of  EMP  by  ultrashort  laser  pulses  in  air  and  on 
dielectric  surfaces  is  analyzed  and  discussed.  Multi-photon  ionization  leads  to  the 
creation  of  a  plasma  filament  in  the  medium.  The  plasma  currents  driven  by  the 
ponderomotive  forces  associated  with  the  laser  pulse  are  the  source  of  the  EMP.  The 
plasma  currents,  however,  decay  rapidly  behind  the  laser  pulse  due  to  the  large  electron 
collision  frequency.  For  a  plasma  density  of  ne  ~  1016  cm-3 ,  the  electron  collision 

frequency  is  ve  ~  5x  1 012  sec"1  (  for  Te  =  1  eV)  and  the  duration  of  the  plasma  current  is 
~  1  /  ve  ~  200fsec .  The  duration  of  the  EMP  at  a  fixed  axial  position  is  roughly  equal  to 
the  laser  pulse  duration  plus  the  plasma  decay  time,  ~  rL  +  1  /  ve ,  which  is  typically  a 

few  hundred  femtoseconds.  The  conversion  efficiency  from  the  laser  pulse  to  EMP 
radiation  is  found  to  be  extremely  small  in  both  air  and  on  dielectric  surfaces.  Typical 
conversion  efficiencies  are  on  the  order  of  1 0“9 .  The  EMP  frequency  spectrum  is  broad 
with  a  characteristic  frequency  determined  by  the  laser  pulse  duration  and  the  electron 
collision  time  and  not  by  the  local  plasma  frequency.  Our  results  are  fundamentally 
different  from  those  in  [30,  31].  For  example,  in  a  steady-state  propagation  regime  the 
EMP  is  non-radiative  and  the  energy  flux  is  axially  directed  with  no  net  radial 
component.  We  show  that  in  a  transient  propagation  regime  radiative  flux  of  EMP 
energy  is  possible  provided  a  Cherenkov-like  condition  is  satisfied.  Several  examples  of 
EMP  generation  have  been  simulated  and  discussed.  Our  results  predict  a  peak  EMP 
power  of  ~  0.8  W  for  the  parameters  of  the  experiment  in  Ref.  17,  corresponding  to  an 
efficiency  of  ~  10'9  for  the  radial  power  flow.  In  another  example,  using  a  high  energy 
laser  pulse  (0.5  J)  we  find  that  the  ionization  losses  of  the  laser  pulse  are  significant  and 
that  the  pulse  eventually  reaches  a  quasi-equilibrium  state.  Finally,  simulations  of  EMP 
generation  from  the  interaction  of  a  laser  pulse  with  a  dielectric  medium  show  that  EMP 
intensities  on  the  order  of  a  MW/cm2  can  be  generated.  Most  of  the  EMP  energy, 
however,  is  reflected  by  the  plasma  formed  at  the  surface  of  the  dielectric. 
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Distorted  Laser  Pulse 


Air 


Figure  1.  Schematic  diagram  for  atmospheric  EMP  generation.  A  ultra-short 
laser  pulse  ionizes  the  air  and  forms  a  plasma  filament.  The  ponderomotive 
forces  associated  with  the  laser  pulse  envelope  drive  radial  and  axial  currents, 
Jr  and  Jz ,  respectively,  which  are  the  source  of  EMP.  The  laser  pulse 
distorts  due  to  ionization  and  the  optical  Kerr  effect.  The  filament  is 
modulated  at  the  plasma  wavelength,  Xp .  The  axial  EMP  electric  field 

outside  the  filament  is  small  due  to  shielding  by  surface  charges  at  the  air- 
plasma  boundary. 
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Figure  2:  (a)  Distribution  of  axial  Poynting  flux  {I z  =  cB2d  /An)  outside  the 
plasma  channel  in  the  (T,r)  plane  for  steady  state  propagation  where  Be  is 
obtained  from  Eq.  (39).  Shaded  contours  are  obtained  for  the  parameters 
neo  =  1-lxlO17  cm'3,  rp  =  0.4  mm  (=  25  dcopo),  rL  =  267  fsec,  XQ  =  0.8  pm 

(Wo/Mpo  =  125  ),  yg  ~  109,  and  IL0  =  7xl013  W/cm2.  Laser  intensity  profile, 

I  id)  i  indicates  that  in  the  steady  state,  the  EMP  duration  is  comparable  with  the 
laser  pulse  duration. 
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Figure  3:  Axially  propagating  EMP  power  (solid  curve),  i.e.,  Pemp(t)  from  Eq. 
(41),  normalized  to  the  peak  laser  power  PL  =  2.8  TW  for  steady  state 
propagation.  Dashed  curve  indicates  longitudinal  laser  intensity  profile  1L  (r) . 
Parameters  are  identical  to  those  of  Fig  2.  In  calculating  PL ,  it  is  assumed  that 

the  laser  spot  size  is  Ra  =  rp  421 ,  where  £  =  8  is  the  multiphoton  ionization  order 
of  Eq.  (5). 
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Figure  5:  Steady  state  EMP  conversion  efficiency  given  by  Eq.  (43)  vs.  laser 
pulse  duration  for  various  laser  pulse  energies. 


Figure  7:  Schematic  diagram  of  intense  laser  beam  incident  on  a  dielectric 
surface. 
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Figure  8:  Comparison  of  longitudinal  electric  field  Ez  obtained  from 
numerical  simulation  (solid  curves)  and  steady  state  analysis  (dashed  curves). 
Plots  show  (a)  normalized  laser  intensity  vs.  T ,  (b)  normalized  electric  field, 
Ez  -  I  <7  I  Ez  /  m  ca)0 ,  vs.  t  at  r  =  0 ,  and  (c)  Ez  vs.  radial  position  at 
t  (opo  /  2n  -  1 .4 .  Laser  parameters  are  Ao  =  0.8  pm ,  rL  =  267  fsec , 

Yg  =13,  and  peak  laser  intensity  IL0  =  7xl013  W/cm2 .  Plasma  parameters 

are  ne0  =  l.lxlO17  cm  3  and  rp  -  0.4  mm .  Simulation  results  are  plotted  at  a 
propagation  distance  of  z  -  6.3  cm  . 
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Figure  9:  Comparison  of  normalized  EMP  magnetic  field  Bg  =  \q\  Bg!mco\ 
obtained  from  numerical  simulation  and  steady  state  analysis.  Panels  (a)  and 
(b)show  Be  in  the  (z,r)  plane  obtained  from  steady  state  analysis  and 
numerical  simulation,  respectively.  Panel  (c)  shows  Be  vs.  normalized  radial 
position  at  z 0)po  /  2/r  =  0.4  obtained  from  steady  state  analysis  (dashed  curve) 
and  simulation  (solid  curve).  Parameters  are  identical  with  those  of  Fig.  8. 
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Figure  10:  EMP  radiance  vs.  normalized  t  at  radial  position  r  -  10  c  /  (Opo 

(outside  the  plasma  filament)  for  the  simulation  of  Fig.  8.  Panel  (a)  shows 
radial  EMP  radiance,  I r  =  cEzB0  /4k  (solid  curve),  and  the  laser  intensity 

profile  I L  (dashed  curve).  Panel  (b)  shows  the  axial  EMP  radiance,  1 2  (solid 
curve),  from  the  simulation  and  the  analytical  calculation  of  I z  =  cB]  14k 
(dashed  curve)  for  steady  state  propagation  where  Bg  is  obtained  from  Eq. 
(39). 
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Figure  11:  (a)  Shaded  contours  of  radial  power  flux,  Ir ,  as  a  function  of 
£  =  z  -  vg  t  and  radial  position,  r ,  showing  radiative  EMP  from  a  modulated 
plasma  filament  with  plasma  density  neo  =l.lxl017  cm-3  and  radius 
rp  =  0.4  mm .  (b)  Ir  vs.  radial  position  at  axial  position  £  =  -  220  pm .  The 
plasma  modulation  is  characterized  by  amplitude  £  =  0.3  and  wavenumber 
K  =  cop  /  Ac  as  defined  by  Eq.  (47).  Laser  parameters  are  X  =  0.8  pm , 

IL0  =  7xl013  W/cm2,  tl  =  267  fsec,  and  yg  =  109 . 
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Figure  12:  (a)  Shaded  contours  of  radial  power  flux,  Ir ,  in  the  (£,r)  plane  at 
a  propagation  distance  z  =  2.2  cm  in  a  modulated  plasma  filament  with 
plasma  density  neo  =  l.lxlO17  cm-3  and  radius  rp  =0.4  mm  .  The  plasma 
modulation  is  characterized  by  amplitude  e  =  0.3  and  wavenumber 
K  =  4  copo  /c  as  defined  by  Eq.  (47).  (b)  Ir  vs.  radial  position  at  axial 
position  £  =  -  220  pm .  Laser  parameters  are  X0  =  0.8  pm , 
ho  =  7xl013  W/cm2,  xL  -  267  fsec,  and  yg  =  109 .  Note  that  the  EMP 

behind  the  laser  pulse  is  confined  within  the  plasma  filament.  The  outward 
radial  flux  observed  within  the  laser  pulse  (-80  pm  <  %  <  0)  is  a 

manifestation  of  transient  interaction  at  the  air-vacuum  interface  at  z  =  0 . 
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Figure  1 3 :  Schematic  diagram  of  the  simulation  of  the  EMP  experiment  of 
Ref.  17.  The  laser  envelope  simulation  (without  EMP)  is  performed  from 
z  =  0  to  z  =  z0 ,  where  z0  =  1 .82  m .  The  full  EMP  simulation  is  carried  out 
for  z  >  zQ. 
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Figure  14:  (a)  Peak  plasma  density  and  laser  intensity  vs.  propagation 
distance  Az  =  z  -  za  for  the  simulation  of  the  experiment  of  Ref.  17.  (b) 

Peak  radial  (dashed  curve)  and  axial  (solid  curve)  Poynting  flux  vs.  Az . 

Initial  (z  =  0)  laser  parameters  for  the  simulation  are  XD  -  0.8  pm , 

PL  =  83  GW ,  zL  =  120  fsec  (FWHM),  and  R0  =  7.5  mm.  Parameters  for  air 
are  taken  to  be  n  -  1  =  2.7xl0'4,  nK  =  3  x  10'19  cm2/W,  U.  =  12  eV, 

^  =  8,  and  I ^  ~1014  W/cm2. 
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Figure  15:  Contours  of  normalized  EMP  fields,  (a)  Ez  =  |  q  \  Ez  /  mca>0 ,  (b) 

=  I  <7 I  Er  /  >  and  (c)  Bg  =  |  g  |  Be  /  wcty0  as  a  function  of 

£  =  z  ~  vg  7  and  radial  position,  r ,  at  propagation  distance  Az  =  10  cm 

corresponding  to  the  simulation  of  Fig.  14.  Unshaded  contours  indicate  the 
laser  intensity  distribution. 
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Figure  16:  Contours  of  (a)  radial  and  (b)  axial  EMP  radiance  as  a  function  of 
=  z  -  v  g  t  and  radial  position,  r ,  at  propagation  distance  Az  =  10  cm 

corresponding  to  the  simulation  of  Fig.  14.  Unshaded  contours  indicate  the 
laser  intensity  distribution. 
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Figure  17:  (a)  Normalized  on-axis  laser  intensity  vs.  time,  (b)  radial  EMP 
Poynting  flux,  lr ,  at  radial  position  r  =  225  pm  vs.  time,  and  (c) 
corresponding  temporal  Fourier  transform  of  Ir  for  propagation  distances 
Az  =  10  cm  (dashed  curve)  and  12  cm  (solid  curve).  The  radial  position 
r  ~  225  pm  is  outside  of  the  plasma  filament.  Simulation  parameters 
correspond  to  those  listed  for  Fig.  14. 
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Figure  18:  Total  EMP  power  emitted  radially  as  a  function  of  time  for  the 
simulation  of  the  experiment  of  Ref.  [17].  Simulation  parameters  correspond 
to  those  listed  for  Fig.  14. 
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Figure  19:  (a)  Normalized  peak  laser  intensity,  I L /IL0  (dashed  curve)  and 
normalized  peak  plasma  density  ne  /  nno  as  a  function  of  propagation  distance 
for  an  initially  collimated  laser  pulse  with  X0  =  0.8  pm ,  PL  =  5.3  TW , 
tL  -  240  fsec ,  and  Ra  =  1  mm.  Parameters  for  air  are  taken  to  be 
no  -  1  =  2.7xl04,  nK  =  3 x  10'19  cm2/W,  U ion  =12eV,  t?  =  8, 

Imp  ~  !0’4  W/cm2 .  (b)  Corresponding  peak  radial  (dashed  curve)  and  axial 
(solid  curve)  Poynting  flux  as  a  function  of  propagation  distance. 
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Figure  20:  Contours  of  laser  intensity  (unshaded)  and  (a)  axial  and  (b)  radial 
Poynting  flux  (shaded  contours)  in  the  (£,r)  plane  at  propagation  distance 
z  =  0.8  cm  corresponding  to  the  simulation  of  Fig.  19. 
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Figure  21 :  Simulation  of  EMP  generation  from  a  dielectric  surface.  Panels 
show  laser  intensity  (unshaded  contours),  and  total  EMP  flux 

I emp  =  dr  +  I2 )u 2  (shaded  contours)  in  the  (z,  r)  plane  at  times  (a)  t  =  0 , 
(b)  t  =  0.3  psec ,  (c)  t  =  0.7  psec ,  and  (d)  t  =  1  psec .  Dashed  line  at  z  =  0 
indicates  the  boundary  between  air  (z  <  0)  and  dielectric  (z  >  0 ).  Laser 
parameters  at  t  =  0  are  -  0.8  pm ,  IL0  =  1013  W I  cm 2 ,  rL  =  200  fsec , 
and  Ra  =0.1  mm.  Parameters  for  the  dielectric  are  taken  to  be  na  =  3 , 
nK=lx  10-15  cm2/W ,  Uion  =  3  eV,  t  =  2,  =  5xl014  W/cm2 . 
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